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AB8TBACT 

This  thetifl  pr«eent8  the  results  of  an  Inrestl^ation  of  the  lateral 
perfonoance  functions  of  an  airplane  by  the  analysis  of  the  transient 
response  of  the  airplane  to  finite  pulse  of  control  surface  defleetioD. 

The  Bethod  is  an  extension  of  pulse  technique  used  by  Clementson, 
Bef.(2),  to  measure  the  longitudinal  performance  functions  of  an  airplane 
and  further  described  by  Seamans ,  Blasingaaa,  ai^  Clementson,  Ref.(3),  for 
more  general  oases.   In  this  thesis  the  successful  extension  to  highly 
oscillatory  lightly  damped  responses  is  shown.  Also  a  possible  technique 
of  analyzing  dirergent  responses  is  indicated. 

The  instrumentation  and  calibration  necessary  to  measure  inputs  of 
aileron  and  rudder  deflection  and  outputs  of  time-rate  of  geometric  yaw  ajjd 
roll  are  described  in  detail.  The  separate  deflections  of  each  aileron 
were  added  electrically  to  gire  the  total  differential  aileron  deflection 
thus  obviating  the  necessity  for  adding  these  deflections  during  analysis. 
Microsyn  signal  generators  were  used  to  measure  aileron  and  rudder  deflections. 
Hate  gyros  were  used  to  measure  the  angular  relocitles  in  roll  and  yaw. 
These  were  the  only  outputs  measured  due  to  the  limited  time  availabl*. 

Determination  of  the  performance  functions  from  the  transient  response 
was  accomplished  by  approximating  the  Toorier  Integrsd  of  the  time  responses 
of  the  inputs  and  outputs  using  the  method  of  summing  approximating  triangles. 

The  technique  of  measuring  the  transient  response  to  finite  pulse 
inputs  is  shown  to  be  accurate,  reliable  and  reproducible  with  a  low  lerel 
of  uncertainties. 

Si>ecial  methods  of  analyzing  lightly  damped  oscillatory  modes  are 
dereloped  and  compared.  The  recommended  method  will  rary  somewhat  wltH  the 
form  of  the  response  curre  but  generally  speaking  the  most  satisfactory 
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method  conalste  of  approilir^tlng  the  Fourier  Traneform  of  the  output  over  the 
early  portion  of  the  time  reaponee  to  a  point  beyond  which  the  response  Is, 
within  the  accuracy  of  raeneurenent ,  only  a  damped  slnueold.   At  this  Talue 
of  time  the  Fourier  Transform  of  the  oscillatory  mode  obtained  analytically 
is  added  vectorially  to  the  Fourier  Transform  of  the  initial  portion  of  the 
response. 

Comparison  with  computed  performance  functions  bpsed  on  calculated 
derivatives  shows  qualitative  agreement  but  clearly  shows  that  these 
computed  performance  functions  may  differ  considerably  from  the  actual 
performance  functions  of  the  full  scale  airplane  in  flight. 

The  saving  in  instrumentation  and  flight  tinie  effected  by  the  pulse 
technique  over  that  required  by  the  sinusoidal  technique  is  a  major 
factor  in  recommending  the  forr.er.  Changes  in  the  mass  and  moment  of 
inertia  due  to  fuel  consurptlon  are  negligible  for  the  short  duration  of 
test  flights  made  possible  by  using  the  pulse  technique. 
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OBJKT 


To  InstnuMnt  a  B-2CJ  airplane  and  neasure  th« 
lateral  transient  reaponsai  in  rate  of  roll  and  rate  of 
^eoaetrio  yav  to  puleee  of  aileron  and  rudder  deflection. 

To  deteniine  the  lateral  perfomanee  fanctlone  of 
the  airplane  froa  thee«  r^eponeee  bj  approxiaat in^  the 
Toarier  Integral!  of  the  input  and  output. 

To  eralnate  this  nethod  of  obtaining  these 
perforaance  functions. 
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IHTRODUCTIOH 

R«cent  deoandt  for  hlgh-performane*  auton»tlc  control  of  aircraft 
for  both  Bllitary  and  coniDerclaJ.  ueee  hare  focu8*d  attention  on  the 
airplane  ae  the  dTnanlc  syetem  being  controlled.   It  is  fundamental  that 
the  eyntheBis  of  a  controlling  aysteii  starts  with  the  beet  poselble 
eraluatlon  of  the  dTnanical  Bysten  to  be  controlled.  Additional  interest 
in  the  dTnaaics  of  the  airplane  arises  from  the  problens  of  stabilizing 
some  of  the  auxiliary  equipment  utilised  in  modem  airplanes,  such  as 
naTigating  amd  fire  control  equipment. 

The  frequency  response  of  an  airplane  can  be  calculated  from  computed 
derlTatlTes  based  on  wind  tunnel  tests.  However ,  the  comparison  of  these 
calculated  responses  with  those  measured  on  the  actual  airplane  has  shown 
appreciable  discrepancies,  due  in  part  to  the  difference  in  stability 
deriratlves  obtained  in  the  wind  tunnel  from  those  of  the  airpljine ,  ajid  in 
part  to  ignorance  of  the  moments  of  inertia  about  the  axes  of  motion. 

The  most  generally  used  methods  of  synthesis  and  design  of  automatic 
control  systems  requires  the  dynamic  characteristics  of  the  airplane  to  be 
presented  as  a  frequency  response.  A  measured  frequency  response  is  most 
often  obtained  by  the  direct  method  of  applying  sinusoidal  inputs  and 
measuring  the  amplitude  and  phase  angle  of  the  outputs.  This  is  done  for 
enough  frequencies  to  define  the  complete  frequency  response. 

Pioneer  work  using  the  sinusoidal  oscillation  technique  was  carried 
out  by  the  Cornell  Aeronautical  Laboratories  with  ontetanding  results. 
These  tests  required  an  auto-pilot,  a  sinusoidal  input  generator  and  measxir" 
Ing  equipment.  Ixtensire  flight  time,  with  its  consequent  expense,  Is 
associated  with  this  technique. 

An  altematlTe  technique,  that  of  measuring  the  trsinsient  response  to 
a  step  function  input,  has  the  drawback  that  such  an  input  to  the  elerators 

-2- 


or  aileroni  will  In  a  Ter7  "iiort  tln«  cAtite  the  RlrplAne  to  rotut*  throng 
a  larg*  angle  and  hence  will  hare  displaced  it  to  far  froB  the  original 
position  of  equili'briuu  that  the  Botion  !•  non-linear.  Such  notion  at  the 
present  tine  can  be  analyzed  only  by  extremely  tedious  and  lengthy  methods 
of  nujnerical  oonputat ions. 

1  technique  which  oVriates  a  large  part  of  the  instrumentation  and  sares 
a  great  deal  of  fll^t  tine  is  that  of  measuring  the  transient  response  to 
a  finite  pulse  of  control  surface  deflection.  This  method  vas  successfully 
used  by  Clementson    *  '   in  determining  frequency  response  in  the  case  of 
longitudinal  motion.  The  mathematicsd  background  and  the  technique  of 

application  are  presented  in  a  paper  by  Seamans,  Blasingame,  and 

-,    ,    (Ref.3.) 
Clementson. 

In  extending  the  pulse  technique  to  the  measurement  of  the  lateral 
■odes  of  motion  of  the  airplane,  the  characteristic  of  spiral  diTergence. 
which  exists  in  many  airplanes,  must  be  taken  into  account.  The  Fourier 
Transformation  method  and  its  approxlaatlon  ®  *  '  are  ralid  only  for 
stable  systems.  Therefore,  it  becomes  necessary  to  identify  and  remor* 
the  unstable  aK>de  prior  to  trans  for  mining  the  pulse  response.  Turthermore, 
the  oscillation  that  occurs  in  lateral  motion  is  lightly  damped.  This 
oscillatory  mode  must  be  giren  special  treatment,  or  the  analysis  will 
become  exceedingly  tedious,  in  the  regions  of  resonance. 

Methods  of  euMdyting  the  oscillatory  mode  were  successfully  used  in 
this  investigation  to  find  the  lateral  performance  functions  of  the  B-25J 
airplane.  However,  this  airplane  did  not  have  any  measurable  characteristic 
of  spiral  divergence,  therefore,  to  illustrate  a  method  of  analysis  when 
such  a  characteristic  is  involved,  a  simple  system  was  analyzed. 
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or  ailerons  will  in  a  rery  short  time  cause  the  airplane  to  rotate  throo^ 
a  large  angle  and  hence  will  hare  displaced  It  so  far  from  the  original 
position  of  eq.alllhrlnB  that  the  notion  Is  noo-llnear.  Sxich  motion  at  the 
present  time  can  be  analyzed  only  "bj   extremely  tedious  and  lengthy  methods 
of  numerical  computations. 

A  technique  which  oVrlates  a  large  part  of  the  Instrumentation  and  sares 
a  great  deal  of  flight  time  Is  that  of  measuring  the  transient  response  to 

a  finite  pulse  of  control  surface  deflection.  This  method  was  successfully 

f Rpf  2  ) 
used  hy  Clementson    *    In  determining  frequency  response  In  the  case  of 

longitudinal  motion.  The  mathematical  "background  and  the  technique  of 

application  are  presented  In  a  paper  hy  Seamans,  Blaslngame,  and 

-.    ^    (Eef.3.) 
Clementson. 

In  extending  the  pulse  technique  to  the  measurement  of  the  lateral 
■odes  of  motion  of  the  airplane t  the  characteristic  of  spiral  dlrergence, 
which  exists  in  many  airplanes*  must  be  taken  Into  account.  The  Fourier 
Transformation  method  and  Its  approxiaation^  e  .  .;  ^^^   ralld  only  for 
■table  systems.  Therefore,  it  becomes  necessary  to  identify  and  remore 
the  unstable  mode  prior  to  transforminlng  the  pulse  response.  Furthermore, 
the  oscillation  that  occurs  in  lateral  motion  is  lightly  damped.  This 
oscillatory  mode  must  be  giren  special  treatment,  or  the  analysis  will 
become  exceedingly  tedious,  in  the  regions  of  resonance. 

Methods  of  analyzing  the  oscillatory  mode  were  successfully  used  in 
this  investigation  to  find  the  lateral  performance  functions  of  the  B-SSJ 
airplane.  However,  this  airplane  did  not  have  any  measurable  characteristic 
of  spiral  divergence,  therefore,  to  Illustrate  a  method  of  analysis  when 
such  a  characteristic  is  involved,  a  simple  system  was  analysed. 


CHAPTXH  1 

DITIHMIHATIOH  07   AIBPLASI  LATIRAL  MOTIOH  PUJOBMAHCI  TUHCTIOKS  BT  PULSE 

TICHHIQUIS 

The  airplane  In  flight  wbj   be  con8ld«r«d.  a*  an  operating  covponent 
which  operate*  on  specified  inputs  to  produce  certain  outputs.  More  thjin 
one  output  aay  result  fron  a  giren  input.  The  operation  on  an  input  to 
produce  a  specified  output  is  expressed  BathBuatically  by  a  concept  of 
wide  generality  called  a  Performance  Operator* 

By  restricting  the  airplane  and  its  controls  to  snail  excursions 
froa  an  equilibrius  condition,  its  motion  may  be  defined  by  a  system  of 
linear  differential  equations.  Then  the  response  can  be  related  to  the 
input  by  means  of  a  Performance  yuaytio^y.  th9   concept  of  Perforaanc* 
Operator  and  Perforaance  Junction  were  both  dereloped  by  Dr.  C.  S.  Draper 
and  his  associates  of  the  Instronentation  lAboratery  of  the  Masesichnsetts 
Institute  of  Technology.  A  complete  discussion  of  their  application  to 
airplane  dynaaics  is  contained  in  Befs.  (l)  and  (lO) .  The  relation  between 
the  two  nay  be  sunaarised  by  the  followini^  equation: 


{PV       .   e''^*  -  e*^^  (pr)r 

'-^in'^out-'  L'lin'Vt- 


Vhen  the  output  of  an  operating  coapoaent  can  be  related  to  the 

input  by  a  set  of  linear  differential  equations,  and  the  input  is 

described  by  a  coaplex  exx>onential,  the  output  is  desoribed  by 

j(w  t  ••■  PA) 
^(out)  ■  ^(out)a  • 

The  ratio  of  the  output  to  the  input  then  defines  the  perfomande 
function  of  the  operating  component  relating  the  two.   In  equation  form 


^"'<  )(vvt' Jt-f»  *  ("'iJ 


«(!»)«  • 


The  tera  a     factors  out,  and  It  is  soon  that  the  performance 
function  is  Independent  of  time .  Tor  uee  In  the  aTnthesls  of  control 
or  ■ta'billilsc  ■jtteBs,  the  perfomanee  fanctlon  is  tLsuallj  presented 
Sraphicalljri  either  as  a  locus  on  the  conplez  plane  or  as  an  anplltude 
ratio  and  phase  auo^le  plotted  rersus  forelni;  an^lar  frequency  on  the 
real  plane.  (l)(2)(3)(l0) 

It  Is  the  need  for  the  frequenej  spectrun  of  the  airplane's  dynamics 
that  has  extended  the  study  of  these  dynaaies  heyond  the  classical  approach 
which,  until  the  last  few  years,  has  been  concerned  with  only  the 
characteristics  equation  of  the  dynaaic  response  emd  sensitlTlties.  This 
approach  does  not  gire  sufficient  information  to  permit  the  desii^  of  a 
hl^--X>erfornance  closed-loop  systea  with  the  airplane  as  the  controlled 
aeaher.  On  the  other  hand,  perforaance  functions  proride  this  necessary 
inforaation  in  a  form  readily  adapted  to  the  prohlems  of  closed-loop  design. 

This  inrestigation  is  concerned  only  with  some  of  the  performance 
functions  which  desori'be  the  lateral  aotiom  of  the  airpleme. 

Specifically,  they  arei 

The  limitation  in  time  arailahle  precluded  the  isrestigation  of  four 
acre  lateral  aotion  perforaance  functions  of  iaportanoe|  e.g. 
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a 

Analjtical  expreBBlona  for  thee*  perfonnBUce  functions  can  "be  o"btalned  from 
the  differential  equations  for  lateral  motion. 

The  Laplace  Traneformatlone  of  the  analytical  equations  are: 

(e-T)v-»-Ur-g^#0 


▼ 

I. 


IZ 

XI      *  *'      "a 

I, 


-tL  >  (-  t~  8  -  L  )r  •►  (•  -  L  )p  =  L.  6  (.) 
T      Iyt      r         p  '^    0   a 


▼I  ♦  (s  -  5^)r  +  (-  —^  B  -  N  )p  =  H.  fi  (b)  -.•  N.  ft  (■)        (l-l) 

Side  force  and  rolling  moment  due  to  rudder  deflection  and  side  forcee 

due  to  rolling  and  yawing  velocities  are  considered  negligible.  The  axes 

and  Bymbols  conform  to  NACA  notation.  Fig.  1-1  pictures  the  azei.  the 

eymholB  are  defined  In  the  Oloeeary,  Appendix  I. 

In  order  to  measure  the  outputs  due  to  a  specific  input,  all  other 

inputs  must  be  zero.  Letting  t     equal  zero  and  specifying  &  as  the 

r  a 

forcing  function,  the  solution  of  equation  (l-l)  yields  the  performance 
function, 

where  V  is  equal  to  r 

This  performance  function  for  the  B-25J  ie  computed  in  Appendix  D,  using 
values  for  the  aerodynamic  derivatives  based  on  wind  tunnel  tests.  Plots 
of  the  amplitude  ratio  and  phase  angle  of  this  calculated  performance 
function  are  given  in  Tigs.  ^-^    and  2-4,  respectively. 

All  of  these  lateral  motion  performance  functions  can  be  obtained 
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•rperimentally  "by  elnuBOldally  oscillating  the  control  aurfao*  and  B«a«'ar- 
log  th«  airplane's  steady  state  responee.  This  sinusoidal  laethod  requires 
a  sine  generator  and  an  auto-pllot  to  supply  the  Inputs,  and  requires  as 
well,  the  usuaI  neaeurlng  equipment.   In  addition,  conslderalile  fll^t  tiae 
Is  needed  to  o'btaln  results  over  a  sufficient  raji^e  of  forcing  frequencies. 

In  the  case  of  nidder  deflection  inputs,  a  step-function  input  would 
yield  a  transient  response  from  vhich  the  perfonn&nce  function  could  be 
detemlned.         *    But  a  step-function  of  aileron  deflection  will 
produce  a  steady-state  rolling  reloclty  which  in  a  short  time  will  generate 
such  a  large  angle  of  .hank  that  the  motion  could  no  longer  be  considered, 
linear. 

The  method  of  using  a  finite  pulse  of  control  deflection  has  the 
advantages  associated  with  obtaining  a  transient  response,  and,  since  the 
control  surface  returns  quickly  to  the  trim  position,  the  disadrantage  of  a 
step-function  input  is  elirelnated.  Finite  pulses  of  elerator  deflection 
hare  been  successfully  used  to  determine  the  airplane's  longitudinal 
performance  functions.        *    But,  to  the  knowledge  of  the  authors, 
this  technique  has  not  been  applied  to  the  lateral  performance  functions  of 
an  airplane. 

The  mathematical  background  for  a  method  of  obtaining  the  performance 

function  from  the  transient  response  to  a  pulse  is  presented  in  Hefs.  (2)  and 

(3).  Briefly,  the  method  consists  of  approximating  the  Fourier  Transforms 

of  the  input  and  the  response.  The  ratio  of  the  latter  to  the  former  is  the 

performance  functlonj  that  is, 

(W)[q^(t)] 
<^'Atq.  .<lo  ]  '   W[q,(t)] 'l-2> 

or,  more  generally. 
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aT)[q^^^(t)] 
^^Uq^.  q^]  "  (LT)[q^^(t)] 


Jw         (1-3) 


When  th«  regponse  Ib  convergent,  It  !«  shown  In  (2)  and  (3)  that 

-Jn^T  w 
(n')[q(t)]'S    (FT)[UT(t)].-   r  q  (t  )•       ^     (I-4) 

n  B  X  |C , . . . 

The  lateral  reeporee  of  an  airplane  requires  special  treatment  'becaus* 
of  the  light ly-daraped  oscillatory  aodc  and  because  of  the  divergent  mode 
found  In  many  airplanes.   The  factored  characteristic  equation  for  the 
lateral  motion  of  the  B-25J,  derived  in  Appendix  D,  shows  both  characteristics: 

[s^  -^  (2)(.14l)(l.34)s  +  (l.34)^][8  +  2.9][(s  -  .0013)]  -  0        (l-5) 

The  lateral  response  of  an  airplane  having  this  characteristic  equation 
would  not  be  convergent.  However,  by  resolving  the  response  into  components, 
the  performance  function  can  be  determined o 

For  example. 


(  >(^ln%ut)       (1-T)[qin(t)] 


> 


The  Laplace  Transform  of  q/   *.Wj4  ^(t)  l-s  obtained  by  first  determining 
the  analytical  expression  for  q/   ^Wjj  \(t).  this  is  obtained  by  plotting 
(on  semilog  coordinates)  the  response,  q/    \(t),  for  large  values  of  time, 
where,  for  practical  purposes,  only  the  divergent  mode  remains.  Trom  this 
plot,  the  slope  and  Intercept  provide  the  necessary  information  to  define 

the  divergent  mode  analytically  and  the  Laplace  Transform  of  this  analytical 

(ref  8) 
expression  can  be  determined.     *     The  Laplace  Transform  of  the  input  it 

equal  to  the  7ourler  Transform,  since  the  input  has  the  value  of  sero  for  all 

values  of  time  less  than  rero  and  its  Integral  is  convergent.   Thus,  the  first 

term  on  the  right  side  of  equation  (I-6)  can  be  found  from  a  time  responee. 

No  evidence  of  a  divergent  mode  was  found  In  this  investigation  for  the 
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B-26J .   Therefore,  It  was  unneceetAry  to  employ  the  technique  described  In 

the  above  paragraph.   Howerer,  to  illuttrate  the  procedure  the  performance 

function  of  an  Idealized  physical  eyetem  with  divergence  was  determined  from 

a  comouted  time  response.   This  performance  function  was  compared  with  the 

ideal  performance  function,  thus  showing  the  feaelbillty  of  the  method.   Thle 

analysis  is  contained  in  Appendix  C. 

The  divergent  component  of  the  tine  response  must  be  subtracted  from  the 

recorded  response.   The  Integral  of  the  remaining  function  of  time  converges. 

However,  btcause  of  the  lightly-damped  oscillatory  mode,  the  approximation 

of  the  Fourier  Transform  by  triangles  becomes  tediously  long.   Thle  problem 

cein  be  overcome  by  subtracting  the  oecillatory  mode.   The  frequency,  damping 

ratio,  and  amplitude  of  the  oscillatory  mode  are  determined  as  described  in 

detail  in  Appendix  D.  This  mode  may  be  either  su"btracted  over  the  region  of 

time  from  zero  to  infinity,  or  it  may  be  subtracted  only  In  the  region  from 

any  selected  time,  t  ,  to  infinity.   In  the  first  case,  the  performance 

function  becomes: 

(yT)[q,    .f        x(t)]-^[q/    V,    Jt  )](rT)[(UlT)(t)]^ 
py  ^        (ontHosc;   -*    ^outMrca^  0      o 

^  ^f^^l'^^  "^  (?T)[UT(t)],  Z  q.,  .(t  ).--^^^1^ 

^  n  =  1,2,3  ^^""^     ° 

-JnAr  w 

The  unit  triangle  of  t  is  here  defined  as  a  triangle  of  unit  height  and  with 

a  base  length  of  2Ar ,   as  in  Ref.  3.   Similarly,  the  unit  right  triaxigle  of 

t  is  defined  as  a  right  triangle  of  unit  height  with  a  base  length  of  At. 

Note  that  a,      ^v/    >.(t)  is  not  the  eame  function  as  q/   ^\/    \(t)  in  eq. 
^(outMrem;  ^(outnrem; 

(1-5). 

In  the  second  of  these  two  methods,  the  performance  function  becomes: 
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(P7) 


^    ^^^(in)*^(out)^ 


P 


(rT)[UT(t)]     L 


n  «:   1 ,2  , . .  .t 


^(out)^^^' 


o 


(1-7) 


(TT)[UT(t)].    I  ^iJ-^J'^'^'^^^" 

n  =  1 ,c ,  t .  • 

If   the  At   aeeoclftted  with  O/      ^x    equals   that   aesociated  with  q.    ,    only  the 

(out;  in 

Buraroatione  remain  in  the  second  tern  on  the  right  eide  of  equation  (l-7). 
However,  if  At^  is  not  equal  to  ^t.,  a  correction  factor  laust  be  applied, 

since  the  Fourier  Transform  of  a  unit  triangle  of  t  is 

(sin  ^jf 
2 


(rT)(UT(t)) 


^    2   ' 


At 


The  Tourler  Tmnsfonns  of  both  a  unit  triangle  and  a  unit  right  triangle,  ae 

functions  of  wAt ,  are  plotted  in  Ref.  (3).  The  Fourier  Transforms  of  both  a 

sine  function  and  a  cosine  function  over  the  region  t   to  infinity  are 

developed  in  Appendix  D,  and  their  application  to  the  determination  of  the 

lateral  motion  performance  functions  measured  in  this  investigation  is 

discussed  In  detail. 

lither  of  the  above  methods  Is  considered  acceptable.  However,  the 

results  of  thin  investigation  Indicate  that  adding  of  the  oscillatory  mode 

at  t   is  a  somewhat  better  method  than  including  the  oscillatory  mode  from 
p  "* 

time  equal  to  rero. 

The  input  pulse  should  have  sufficient  area  to  produce  a  desirable 
magnitude  of  response,  but  should  not  cause  motion  beyond  the  linear  region 
of  the  airplane's  aerodynamic  characteristics.   Also,  to  simplify  the  data 
reduction,  it  is  desirable  that  the  pulse  be  of  short  duration.   The  pulses 
used  in  this  investlp^ation  were  roughly  estimated  from  previous  data  on 
static  sensitivities,  such  as  p  b/2U,  Ref.  4,  together  with  the  experience 
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gained  from  the  longitudinal  pulse  teeti  on  thi^  airplane.  The  llnearltj 
ranges  of  the  litxkagee  for  the  control  deflection  plcfc-offs  (mlcroeyn  signal 
generators)  and  the  sensltlrltles  of  the  rate  gyros  were  detemlned  on  the 
basis  of  those  estinates  with  a  100  percent  safety  factor.  After  flight 
trials,  the  rough  estimate  of  the  aileron  pulse  was  increased  by  about 
50  percent  in  oagnitude,  with  the  duration  remaining  as  originally 
estimated.  On  the  other  hand,  the  estimated  rudder  pulse  was  reduced 
slightly  after  flight  trials. 
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CRAPTIR  2 
IHSTBUMXHTATIOK  USD   FLIGHT  TIST  PHOCIDUHI 
A.   INSTHDMKHTATIOH 

The  airplane  employed  ab  a  te«t  rehlcle  for  thi»  investigation  was  a 
USAT  North  American  B-25J ,  Army  Number  44-30328.   The  airplane  had  been 
stripped  of  all  araament .  A  photograph,  a  three-view  drawing  and  a  table 
of  general  apeclf icatione  and  dimensione  are  presented  in  Appendix  A.  This 
airplane  had  previously  been  instrumented  to  record  the  normal  acceleration 
and  pitching  velocity  responses  to  a  recorded  pulse  of  elevator  dlsplacempnt . 
The  amplifying  and  recording  systems  used  in  the  earlier  testing  were  retained 
for  the  present  studies  augmented  by  the  following: 

1.  Adjustable  rudder  stop  device  (mounted  on  the  co-pilot's  right 
rudder  pedal),  fulfilling  the  dual  purpose  of  providing  rudder 
pulses  of  any  desired  fixed  magnitude,  and  of  locking  the  rudder 
while  aileron  pulses  were  being  applied.   See  Figs.  4  and  5, 
Appendix  B . 

2.  Adjustable  aileron  stop  device  (mounted  on  the  co-pilot's  control 
column) ,  fulfilling  the  dual  purpose  of  providing  aileron  pulses 
of  any  desired  fixed  magnltud*.  and  of  locking  the  ailerons  while 
rudder  pulses  were  being  applied.   See  Figs.  1,  2  and  3,  Appendix  B. 

3.  Mlcrosyn  signal  generator  unite,  linked  to  each  aileron  for  recording 
the  magnitude  of  the  aileron  displacement.   See  Figs.  6  and  7, 
Appendix  B . 

4.  Mlcrosyn  signal  generator  unit  mounted  at  the  inidder  control  bell 
crank  (right  rudder)  for  recording  the  magnitude  of  the  rudder 
displacement.  See  Fig.  8,  Appendix  B, 

5.  Eate  gyro,  for  measuring  the  angular  velocity  in  yaw.   See  Fig.  15, 
Appendix  B  . 

6.  Hate  gyro,  for  measuring  the  angular  velocity  in  roll.   See  Fig.  15, 
Appendix  B . 

7.  llectronic  current  regulator,  used  to  regulate  a  set  value  of  stiffness 
current  applied  to  the  elastic  restraint  generators  of  the  roll  and  yaw 
gyros. 

Four  channels  of  amplification  and  four  channels  of  a  Consolidated 

recording  oscillograph  were  utllited  to  record: 
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1.  Aileron  Input  pulae*. 

2.  Rudder  Input  pulsus. 

3.  Roll  angular  velocity. 

4.  Taw  angular  velocity. 

The  ampllfiert  and  oscillograph  are  described  In  Appendix  B. 

Several  ground  callbratione  were  made  before  any  flights  were 
attempted,  and  a  final  ground  calibration  was  performed  prior  to  the 
final  flight.  The  methods  of  calibration  and  the  calibration  curves  are 
presented  in  Appendix  B. 

The  center  of  gravity  of  the  airplane  was  determined  by  means  of  a 
USAT  electronic  weighing  kit.  The  center  of  gravity  of  the  airplane  was 
uaintained  near  the  center  of  gravity  of  the  front  main  fuel  tanks  by 
shutting  off  the  other  tanks  during  flight.  This  caused  the  shift  In  the 
center  of  gravity  during  flight  to  be  negligible. 

An  accurate  determination  of  the  moments  of  inertia  of  the  B-25J 
airplane  for  use  in  the  theoretical  calculations  of  this  thesis  posed  a 
difficult  problem.  Obviously,  they  could  be  determined  most  accurately 
by  experimental  means,  but,  because  of  the  site  and  weight  of  the  B-26J, 
such  a  procedure  was  beyond  the  scope  of  the  present  study.  As  a  practical 
alternative,  the  moments  of  inertia  were  fixed  by  careful  estimates  from 
known  data  Already  established  on  similar  aircraft. 

Estimates  were  first  based  on  the  fairly  complete  data  for  a  B-25J 
contained  in  Ref.  1.  However,  a  subsequent  comparison  of  the  values  for 
lyy  contained  in  the  data  with  the  value  for  I__,  determined  for  our  test 
airplane  by  other  investigators  (Ref.  2)  Indicated  that  the  test  airplane 
and  the  reference  airplane  differed  somewhat  in  configuration.   Therefore, 
the  final  estimates  were  based  on  a  comparison  with  similar  data  from  the 
A-26  airplane  discu^ned  in  Ref e .  7  and  9.   This  comparison  was  expected  to 
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Increase  the  validity  of  the  eatlirated  figures  because  the  test  airplane 
and  the  A-26  had  similar  dimensions  and  differed  only  slightly  In  weight. 
This  comparison  yielded  somewhat  higher  figures  for  I__  and  I^^  than  those 
given  in  Ref.  1,  in  line  with  the  higher  value  for  I_Y<letermlned  experiment" 
ally  for  the  test  airplane,  to  which  reference  has  already  been  made, 
B.   FLIGHT  TEST  MITHODS 

Thirty  minutes  prior  to  a  flight,  all  equipment  was  tum»d  on 
temporarily  and  the  following  items  were  checked  to  Insure  proper 
operation: 

1 .  300-Tolt  d-c  supply 

2.  26-Tolt  three-phase  gyro  wheel  excitation 

3.  150-volt  d-c  supply 

4.  Oyro  stiffness  current 

5.  Amplifier  outputs 

6.  Paper  supply  in  oscillograph 

7.  Frequency  of  the  Inverter  output 

Soon  after  take-off,  all  equipment  was  turned  on  a^ain  —  this  time 
for  the  duration  of  the  flight  —  to  insure  sufficient  warn  up  time,  and 
the  Items  mentioned  above  were  rechecked.  Then,  with  the  airplane  trinmed 
at  175  mph  emd  10,000  ft.  altitude,  test  pulses  were  applied  so  that  the 
sero  position  of  the  light  traces  of  the  oscillograph,  the  amplifier 
attenuation  settings,  and  the  input  pulse  amplitude  could  be  adjusted, 
▲n  effort  was  made  to  keep  the  attenxiation  settings  constant  so  as  to 
reduce  the  number  of  calibrations  required. 

After  the  preliminary  adjustments  had  been  made,  the  airplane  was 
:5  carefully  trimmed,  with  all  controls  held  firmly  against  the  stops,  the 
pilot  then  started  the  oscillograph,  and  three  seconds  later  gave  the 
signal  for  thp  Impulse.   The  oscillograph  was  stopped  by  the  oscillograph 
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operator  when  the  oscillatlone  had  attenuated  to  a  negllgl'ble  aniplltude. 

The  number  of  control  stops  used  In  these  tests  made  It  necessary  for 
the  pilot  and  co-pllot  to  follow  a  certain  sequence  of  operations.   This 
sequence  was: 

1.  Prior  to  recorded  run: 

a.  Pilot  trinmed  airplane  to  maintain  straight  cuad  level  flight, 
holding  the  elevator  control  against  the  elevator  stop. 

"b.  Co— pilot  adjusted  aileron  and  rudder  stops  to  give  desired  pulse 
amplitude*.  (This  was  done  by  applying  pulses  of  various 
magnitudes,  within  the  linear  range,  until  the  desired  response 
was  observed  on  the  oscillograph.) 

2.  For  a  recorded  run: 

a.  Pilot  triiimed  airplane  in  straight  and  level  flight,  holding  the 
,  elevatoi*  control  against  the  elevator  stop. 

b.  Co-pilot  moved  adjustable  trim  screw  for  rudder  until  the  proper 
stop  was  engaged  I  the  rudder  was  then  held  firmly  against  this 
stop. 

c.  Co-pilot  moved  the  adjustable  trim  stop  pin  for  the  aileron  until 
the  proper  stop  was  engaged,  the  aileron  control  was  then  held 
firmly  against  this  stop. 

d.  Pilot  called  oscillograph  operator  and,  if  the  recording  equip- 
ment was  ready,  pushed  the  button  to  start  the  oscillograph. 

•.  Three  seconds  after  starting  the  oscillograph,  the  pilot 
signalled  the  co-pilot  to  apply  the  pulse. 

On  each  flight,  two  records  were  made  for  each  of  the  following  types 

of  pulse: 

1.  Positive  aileron  (right  aileron  down) 

2.  Negative  aileron 

3.  Right  rudder 

4.  Left  rudder 

After  the  above  records  had  been  taken,  the  airplane  was  held  on  a 
straight  course  at  10,000  ft.  for  a  sufficient  time  to  make  calibration 
records.   The  method  of  calibration  is  discussed  in  detail  in  Appendix  B. 
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In  addition  to  the  oscillo/^rftph  records,  the  following  data  were  also 
recorded  for  each  run: 

1.  Attenuation  settings 

2.  Direction  of  control  deflection 

3.  Altitude 

4.  Indicated  airspeed 

5.  Outside  air  temperature 

6 .  Fuel  load 

7.  Inline  rpm 

8.  Inverter  frequency 

9.  Gyro  elastic  restraint  current 

10,  Signal  generator  excitation  current 
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CHAPTIB  3 
BISULTS  AND  DISCUSSION 
A  satiefactory  Instrumentation  and  calibration  of  the  test  vehicle 
wa«  achieved.  It  is  shown  in  Appendix  B  that  the  maiiBtun  oncertainty  in 
calibration  was  1.7  percent. 

The  flight  test  data  was  reduced  by  the  methods  outlined  in  Eef .  3 
with  the  following  results: 

(a)  The  uncertainty  in  reproducing  results,  obtained  by  analyzing 
three  independent  records  of  rolling  velocity  response  to  input  aileron 
pulses,  was  fouac!  to  be  1.5  percent  average  deviation  from  the  average 
measured  performance  function  in  amplitude  ratio,  and  2  average  deviation 
in  phase  angle,   (See  Figures  3-1  and  3-2), 

(b)  The  comparison  of  the  experimentally  determined  performance 
function  with  the  performance  function  calculated  from  wind  tunnel  stability 
derivatives  for  rolling  velocity  response  to  aileron  Input  showed  very  close 
eigreement  in  phase  angle  at  all  frequencies  within  the  range  Investigated 
(0,5  to  6.0  rad/sec).  The  comparison  in  AH  showed  both  curves  following 
the  same  general  pattern  with  the  calculated  performance  function  having 

a  lower  AH  at  all  frequencies  except  the  resonant  frequency.  The  resonant 
frequenclti  were  displaced  both  in  amplitude  ratio  and  In  frequency  In 
accordance  with  the  respective  differences  in  damping  ratio  and  undamped 
natural  frequency  of  the  oscillatory  modes.   (See  Figures  3-4  and  3-2). 

(c)  The  eensltivities  of  the  calculated  and  experimental  performance 
functions  for  rolling  velocity  response  to  aileron  input  compared  as 
follows: 

a  A. 

8(A) (d  W  ),     w      --1-00  1/eec 
a  X  (exper Haver) 
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2U 
addltioHAl  eenfiltiTlty  Is  available,  I.e., 


Utilizing  ^  vs  8   data  obtained  from  flight  teet  data  (Ref .  4)  an 
cU      a 


a  X  (fit .  tcBt  data; 
(This  sensitivity  allows  a  rough  check  on  overall  calibration  of  the 
InBtrunentation  eyBtem) . 

(d)  Comparieon  of  the  experimental  performance  function  for  yawing 
rate  out  to  aileron  input  with  the  calculated  performance  function  shows 
good  correlation  in  AH.  The  PA  comparison  is  satisfactory  for  low 
frequencies  (.5  to  3  rads/eec) .  However,  At  high  frequency  a  wide  scattering 
of  points  is  produced  using  three  different  systems  of  analysis.   (Ses 
Tigures  3-3  and  S-^). 

(e)  Sensitivity  comparison  of  the  calculated  and  experimental  perform- 
ance functions  for  yawing  velocity  response  to  aileron  input  is  as  follows: 

(f)  From  an  investigation  of  three  different  methods  of  dividing  the 

output  into  comr/onents,  as  discussed  in  Appendix  D,  it  was  found  that  the 

most  accurate  and  fastest  system  of  graphical  analysis  for  use  in  the  case 

of  a  lightly  damped  oscillatory  response  was  method  (2)  of  Table  2, 

Appendix  D.   In  this  method  the  usual  triangular  approximation  is  used 

until  the  time  response  is  only  a  pure  damped  sinusoid.   At  this  point, 

t  =  t  ,  the  Fourier  Transform  of  a  cosine  or  sine,  whichever  is  appropriate, 
P 

is  added. 

(g)  Investigation  of  the  scattering  of  points  at  high  frequency  in 

the  performance  function,  (^)/4W*  y  )'    showed  that  the  uncertainty  was 

a  w 
a  function  of  the  following: 
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(1)  Th«  m«R*iareinf>nt  uncertainty  in  detemlning  rector  magnitudes  from 
th«  Oflcillograph  record, 

(2)  Deteraination  of  the  proper  bate  line. 

(3)  The  time  response  pattern  —  depends  on  whether  the  aaiimtUB  ordinate 
occurs  early  or  late  in  the  tlA*  response,  and  on  how  lar^e  a  proportion  of 
the  response  is  purely  oscillatory. 

(4)  The  ratio  of  magnitudes  of  the  last  plotted  rector  in  the  output 
BUHimation  to  the  resultant  rector  representing  the  response. 

(5)  The  method  of  component  "breakdown. 

(6)  The  frequency  —  uncertainty  increases  with  frequency  in  the  usual 
case. 

(7}   Size  of  At  used  in  the  approximation  of  the  input  and  output. 

Inetrunentation  and  the  actual  test  flights  for  this  project  consumed  a 
relatirely  large  percentage  of  the  time  arailahle .  Consequently,  the  time  re- 
maining for  analysis  of  the  oscillograph  records  required  that  only  those 
analyses  most  pertinent  to  the  objectires  of  the  project  be  undertaken.  Since 
the  main  objectire  of  the  project  was  to  inrestigate  the  practical  aspects  of 
the  application  of  the  pulse  technique  to  dynamic  systems  in  general,  and  since 
determination  of  the  performance  functions  of  the  test  rehicle  were  of  secondary 
importance,  the  scope  of  data  analysis  was  limited  to  the  following  inreetigation: 

1.  Determination  of  the  uncertainty  in  repeating  results  by  obtaining 

a  performance  function  from  each  of  three  oscillograph  records  independently, 
using  aileron  pulse  Inputs  and  rolling  relocity  response. 

2.  Ivaluation  of  sereral  different  graphical  plotting  methods  as  to 
accurftcy,  speed,  and  applicability  to  various  response  configurations. 

3.  Conqparisoh  of  the  experimental  performance  function  with  the 
calculated  performance  functions,  using  response  in  rolling  relocity  and 
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7«wlng  velocity  to  aileron  inputs  as  exeunples. 

4.   Seeking  out  of  liBltatlons  Impoeod  on  the  accuracy  of  the  pulse 
technique  due  to  the  approximations  involved  In  the  process  of  converting 
froB  the  time  domain  to  the  frequency  domain. 

Table  1  of  Appendix  D  lists  the  oscillograph  records  analysed,  and  the 
method  of  analysis  used. 

The  uncertainty  in  reproducing  results  using  the  pulse  technique  is 
considered  to  he  entirely  satisfactory,  and  well  within  the  demands  of  usual 
engineering  requirenents .   In  evaluation  of  the  ability  of  the  pulse 
technique  to  reproduce  results  two  investigations  were  pertinent.  They 
were:   (l)   the  pbility  to  produce  the  eane  performance  function  by  analyzing 
different  records,  and  (2)   the  ability  of  two  different  analysis  methods  to 
produce  the  same  performance  function  when  applied  to  the  same  oscillograph 
record.  The  results  of  the  first  investigation  are  shown  in  Figures  3-1  and 
8^2.  The  three  oscillograph  records  used  to  obtain  the  experimental  results 
shown  in  Figures  5-1  and  3-2  are  tabulated  in  Table  1,  Appendix  D.  The 
experimental  curve  shown  Is  the  average  of  three  performance  functions  deter- 
mined independently.   The  maximxim  AH  deviation  from  the  mean  curve  is  5.9 
percent  at  u^^  =  1.5  rad/sec.,  and  the  average  AR  deviation  from  the  mean 
curve  is  1,5  percent.  For  the  phase  angle,  the  maxinum  deviation  is  6 
and  the  average  deviation  is  2  . 

The  uncertainty  in  repeating  results  by  utilizing  the  three  graphical 
methods  listed  in  Table  2  of  Appendix  D,  was  found  to  be  a  variable  depending 
on  the  response  pattern  in  the  time  domain,  and  the  forcing  frequency.   In 
general,  the  method  described  as  Method  (2)  in  Table  2,  Appendix  D,  offers 
the  least  uncertainty.   This  method  utilizes  the  same  degree  of  approximation 
in  determining  the  Fourier  Transform  of  the  input  and  output  from  the  records 
in  the  time  domain  out  to  t  =  t  .  Disregarding  frequency  effecti,  the 

-21- 


uncertainty  of  this  method  nay  be  held  fairly  constant  with  Tarlatlon  in 

reeponae  pattern  by  selecting  t  at  increasingly  higher  Talues  as  the 

response  pattern  becomes  increasingly  more  oscillatory.  Comparing  the 

pattern  of  rolling  relocity  response  with  the  pattern  of  yawing  relocity 

response  using  aileron  inputs  for  both,  it  was  found  that  in  the  former 

a  low  proportion  of  the  response  was  oscilletory  while  in  the  latter  a  high 

proportion  of  the  response  was  oscillstory.  Consequently,  in  the  latter 

case,  if  t  was  not  increased  beyond  the  t  used  in  the  rolling  response 
P  P 

case,  a  higher  proportion  of  the  response  was  represented  by  a  Tourier 

Transfora  which  is  exact,  while  the  input  remains  an  approximation  to  the 

exact  Tourier  Transform.  The  ultimate  limit  is  to  let  t  increase  to  a 

P 

time  beyond  which  the  response  la  negligible,  plotting  all  yectore  directly 

from  the  record.  This  is  listed  in  Table  2,  Appendix  D  as  method  (l) . 

Thie  method  offers  the  least  error  in  the  hypothetical  case  where 

instrumentation  inaccuracies  and  noise  are  at  such  a  low  lerel  that 

ordlnates  may  be  measured  accurately  at  large  ralues  of  time. 

Increasing  t  as  the  percentage  of  oscillatory  portion  of  the  responee 
P 

increases  is  also  important  in  nullifying  errors  occassioned  by  choosing  t 

P 

at  a  point  which  is  not  exactly  a  true  peaic  or  lero  of  the  response  ctirre. 

An  error  of  this  nature  eerret  to  shift  the  phase  angle  of  the  added  vector 

for  the  oscillatory  remainder  (t  =  t  )  by  an  amount  equal  to  uj^t  .   If  the 

P  f 

magnitude  of  the  last  rector  plotted  in  the  summation  of  triangles  is  only 
a  small  percentage  of  the  added  rector,  the  error  in  phase  angle  may  be 
large.  However,  if  the  magnitude  of  the  vector  added  is  of  exactly  the  same 
length  as  the  last  rector  of  the  summation,  an  error  in  t  as  large  as  t  ■  At 
will  yield  no  error  in  the  PA  of  the  response  rector. 

The  method  of  dividing  the  output  into  components,  described  in  Appendix 
D,  Table  2,  as  method  (3),  appears  to  be  subject  to  the  greatest  inaccuracies 
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any  of  the  methfxlt  InreBtl^ated  In  analyzing  a  lightly  damped  oecillatory 
reeponse.   Theoretically,  this  method  of  analyeia  le  ae  sound  as  any  of  the 
methods  investigated.   Howerer,  the  practical  application  of  this  method 
preeentB  an  additional  source  of  error  since  the  Input  and  output  are  not 
represented  with  the  same  degree  of  approximation.   In  the  lightly  damped, 
highly  oscillatory  response,  essentially  the  entire  output  is  represented 
by  the  Fourier  Transform  of  the  oscillatory  mode,  which  is  the  exact 
transform  of  an  analytical  expression.   On  the  other  hand,  the  input  is 
approximated  "by  a  vector  addition^  the  exactness  of  which  depends  on  the 
accuracy  of  curve  fitting  achieved.   At  relatively  high  frequency,  the 
Bmall  resultant  vector  of  the  remainder  summation,  Although  approximated 
with  the  same  degree  of  exactness  as  the  input  resultant  vector.  Is 
Ineffectual  in  changing  the  total  output  phase  angle  which  Is  dominated  by 
the  large  vector  representing  the  oscillatory  mode.   Consequently,  the  output 
phase  angle  becomes  nearly  static  as  the  limiting  phase  angle  of  the  oscillatory 
mode  is  approached,  while  the  input  phase  angle  continues  to  grow  as  the 
forcing  frequency  increases. 

In  performing  the  component  breakdown  necessary  in  this  system,  additional 
measurement  errors  may  accumulate  if  more  than  one  set  of  measurements  are 
taken  from  the  response.  That  is,  in  removing  the  oscillatory  mode,  the 
nrdlnates  as  determined  by  the  analytic  expression  of  the  oscillatory  mode 
should  be  subtracted  algebraically  rather  than  graphically  from  measured 
ordlnates  of  the  total  response.  lurthermore,  since  the  oscillatory  mode 
represents  the  major  part  of  the  response,  great  care  must  be  exercised  in 
determining  W  ,  DR,  and  <\,    v. 

In  Figures  3-3  and  3-4,  the  results  of  analyses  by  the  various  methods 

discussed  are  shown.   The  scatter  of  points  at  high  frequency  shows  the 

inaccuracies  of  the  several  methods.   Theoretically,  the  PA  curve  should 

level  off  in  the  vicinity  of  -90  phase  angle.   However,  Fig.  3-4  shows  the 
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•xperlraentAllj  determined  PA  lagging  leee  at  ue.  of  5  and  6  rad/sec  than 

ftt  3  rad/eec.   In  Inr^ttigatlng  this  diecrepancy  between  theoretical  and 

experimental  results ,  the  question  of  accuracy  in  approximating  the  Fourier 

Transform  arose. 

To  determine  how  closely  the  approximation  of  input  and  output  time 

response  "by  summation  of  triangles  fits  the  actual  curree,  frequency  spectra 

of  yaw  response  and  aileron  input  were  prepared  by  numerical  means.  The 

results  are  shown  In  Figures  3-5  and  3-6.  The  circled  points  show  the 

results  of  the  numerical  analysis  and  the  dotted. curre  shows  the  mean 

curre  of  the  experimental  points  obtained  in  Figures  3-3  and  3-4.  The  phase 

angles  determined  by  the  numerical  analysis  show  a  variation  of  more  than  30 

from  the  arerage  curre  at  many  points.  The  explanation  for  such  a  discrepancy 

was  found  to  lie  in  a  difference  in  measureBwnt  of  the  ordinates,  [q(t  )], 

between  those  used  in  the  numerical  analysis  and  those  used  in  the  graphical 

ainalysit.  The  measurements  of  the  reipectire  areas,  for  each  analysis, are 

as  follow: 

ARIA  AB£A 

Numerical  Graphical 

Analysis  Analysis 

1   positive  lobe  0.268  in  see       0.26  in  sec 

St 

1  negative  lobe  0.6170  in  sec      0.5990  in  sec 

nd 

2  positive  lobe  0.1126  In  sec      0.1096  In  sec 

In  the  two  most  important  lobes,  the  deviation  in  measurement  of  area 
is  approximately  3  percent,  emd  in  each  case  the  graphical  analysis  is 
the  smaller  area.  Thus  it  is  noted  in  Figure  3-7  that  a  small  uncertainty 
of  3  percent  in  area  determination  may  cause  as  much  as  30  difference  In 
phase  angle,  provided,  the  Tector  pIo.^  of  the  output  '  'winds  up* '  near  the 
origin.  Since  the  numerical  analysis  gives  only  a  cheek  on  how  closely 
the  triangular  summation  approximates  the  area  of  the  curve,  the  performance 
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function  wftfl  a^aln  determined  At  Bereral  forcing  frequenclta  UBln^  the  fl&me 
ordinatea  as  used  In  the  nuaerlcal  ajial78ls.   The  trianglee  in  Tigore  2n6 
•how  the  points  thus  determined.   It  is  seen  that  the  uncertainty  of  repro~ 
ducin^  results  is  satisfactory,  indicating  that  the  use  of  triangles  in 
approximating  the  area  yields  approximately  the  same  accuracy  as  the  nuaerlcal 
approximation.   It  is  now  apparent  that  any  insiccuracies  in  seasuring  the 
ordinate*  h«re  %  critical  effect  on  the  performance  function  when  the  end 
point  of  the  vector  summation  is  near  the  origin.  However,  the  PA  variations 
due  to  inaccuracies  in  measurement  of  ordinates  is  not  constant,  hut  is  a 
function  of  frequency.  Although  the  measuring  error  in  this  case  gives  an 
increase  of  ordinate  magnitude  in  both  the  positive  and  negative  lohes,  the 
errors  in  a  positive  lohe  do  not  nullify  the  errors  in  a  negative  lobe 
except  at  certain  frequencies  which  are  indeterminate  unless  the  actual 
measurement  errors  are  known.  Figure  3r>7  shows  the  results  of  measurement 
error  at  u).  ■  5  rad/sec.  Due  to  the  phase  angle  shift  with  each  added 
vector  (w^t)  the  errors  are  apportioned  in  accordance  with  the  sine  and 
cosine  of  the  vector  phase  angle  multiplied  by  the  magnitude  of  the  error 
in  that  vector.   If  the  sum  of  the  imaginaries  and  reads  of  the  measurement 
errors  add  up  to  give  a  vector  which  has  an  angle  corresponding  to  that  of  the 
original  resultant  vector,  no  change  in  PA  will  result.  However,  this  rarely 
happens,  and  at  high  frequencies  a  large  phase  angle  change  results.  Hote  in 
Tig.  3-7  that  the  largest  differences  in  the  vector*  of  the  first  two  lobes 

0 

occur  in  measuring  vectors  ^ich  have  phase  angles  close  to  ■♦^0  ,  i.e.,  vector 
10  in  the  first  lobe  and  vectors  28  and  30  in  the  second  lobe .  All  of  these 
difference  tend  to  increase  the  phase  angle  at  »-  ■  5.  lach  frequency  will 
have  a  different  total  PA  error  due  to  measurement  uncertainties,  since  the 
phase  angle  of  each  vector  changes  with  each  different  forcing  frequency. 
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Again  not*  th&t  measurement  error*  hare  an  intignif leant  effect  unleea 
the  end  point  of  the  BuiinnRtlon  is  near  the  origin.  At  »-  =  1.76  rad/sec 
the  PA  is  the  same,  whether  determined  graphically  or  numerically,  and  the 
effect  of  loeasurement  differences  on  uncertainty  in  determining  the 
performance  function  is  negligible  as  shovn  by  the  superiaaposed  point* 
lying  on  the  arerage  experimental  curve. 

In  general,  it  appears  that  any  uncertainty,  whether  in  measurement 
of  ordinates  or  base  line  location,  gires  a  magnified  uncertainty  in  PA 
if  the  end  point  of  the  suaoDBtlon  is  near  the  origin. 

The  phase  angle  uncertainty  is  not  due  to  a  breakdown  in  the  theory  of 
the  pulse  technique,  nor  caused  directly  by  the  triangular  approximation 
of  the  Input  and  response.  Bather  it  is  a  function  of  uncertainty  in 
measurement  of  ordinate*  and  in  plotting  which  accumilate  ui^er  the 
condition*  of  high  frequency  (rapid  phaee  angle  change  of  each  added 
▼ector)  and  low  amplitude  ratio  (require*  and  point  to  be  near  the  origin). 
TtLrthermore ,  this  uncertainty  in  phaee  angle  is  not  restricted  to  the 
predominately  oecillatory  re*pon*e.  It  merely  appears  at  a  lower  frequency 
in  this  case  due  to  the  rapid  attenuation  in  amplitade  ratio  (small 
Tector*  in  the  early  portion  of  the  respon**) .  The  ease  phaee  angle 
uncertainty  would  occur  in  the  rolling  relocity  reponee  to  aileron 
input  if  the  perfomence  function  detemlnAtion  had  been  carried  out  to 
thoee  high  frequencies  whore  the  lunplltude  ratio  wa*  greatly  attenuated* 
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CHAPTKE  4 
CONCLUSIONS  AND  HECOMMKNDATIONS 
A.   CONCLUSIONS 

As  a  result  of  this  investigation  the  following  conclusions  are 
submitted: 

1.  The  pulse  technique  is  a  practical  and  suitably  accurate  method 
for  determining  aircraft  performance  functions  from  flight  tests. 

2.  Iconony  in  flight  time  and  instrumentation  are  gained  by  the 
pulse  technique  as  compared  to  the  pinusoidal  method. 

3.  The  accuracy  of  the  results  obtained  is  enhanced  by  the  short 
duration  of  test  flights  in  that  the  changes  in  mass  and  mooents  of 
inertia  due  to  fuel  consttmptlon  are  negligible, 

4.  Irroneous  performance  function  phase  angles  will  in  general  be 
obtained  at  relatively  high  frequencies  when  the  amplitude  ratio  has 
attenuated  to  a  very  small  value  (order  of  five  percent).  This  is  due  to 
the  fact  that  the  errors  in  measurements  of  the  ordinates  of  the  recorded 
response  accumulate  at  those  frequencies  where  the  series  of  vectors, 
representing  the  ordinates  and  their  individual  phase  angles,  encircle 
the  origin, 

5.  Lightly  damped  oscillatory  responses  may  be  analyzed  with  reasonable 
simplicity  and  accuracy.   However,  in  the  case  of  a  response  which  is  almost 
wholly  oscillatory,  uncertainty  in  the  data  may  cause  anomalous  phase  effects 
to  appear  in  the  plots  of  the  performance  function  at  the  higher  forcing 
frequencies  when  the  amplitude  ratio  is  considerably  attenuated.   These  phase 
effects  are  more  pronounced  in  the  performance  functions  which  are  highly 
dependent  on  coupling  terms  in  the  equations  of  motion  rather  than  the 
forcing  function  Itself,  since  high  attenuation  in  AR  occurs  at  a  relatively 

low  frequency.   For  example,  PFr  ^  ^  ^  and  PFj-     -.  are  such  performance 

a^         rX 
functions. 

-34- 


6.  A  predominately  oeclllfttory  response  le  best  approximated  by  use 

of  iBosceles  triangles  to  a  value  of  time  beyond  which  only  a  damped  sinusoid 
reBAlns.   The  Fourier  Transform  of  the  sinusoid  is  then  obtained  by  trans- 
forming the  analytic  expression  which  is  then  added  vectorially  to  the 
approximate  transform  of  the  early  portion  of  the  response. 

7.  The  most  important  sources  of  error  in  reduction  of  data  pre: 

a.   Choosing  the  base  line  on  the  oscillograph  record. 
■fa.  Measurement  of  the  ordlnates  from  this  base  line. 

8.  The  calculated  and  experimental  performance  functions  are 
qualitatively  comparable,  but  the  quantitative  differences  indicate  the 
requirement  that  the  dynamic  charpcterietics  of  all  high  speed  and  non- 
conventional  aircraft  be  determined  by  flight  tests. 

B.   EKCOMMENDATIONS 

To  aid  further  research  in  the  use  of  the  pulse  technique  to  determine 
aircraft  performance  functions  the  following  recommendations  are  presented: 

1.  Obtain  relatively  long  records  of  lightly  damped  oscillatory 
responses  to  facilitate  the  determination  of  W  and  DR  of  the  oscillation. 

2.  Amplify  low  amplitude  predominately  oscillatory  responses  as  tnach 
as  is  consistent  with  the  noise  level  and  linearity  of  the  recording 
instruments,  the  friction  torque  level  of  the  gyros,  and  the  linearity  of 
aircraft  motion. 

3.  Limit  the  shock  mounting  of  the  gyros  to  landing  and  take-off 
only.   During  tests  the  gyros  should  be  rigidly  attached  to  the  aircraft 
structure  so  that  the  input  to  the  gyros  is  the  true  angular  velocity  of 
the  Airplane  and  is  not  modified  by  the  flexibility  of  the  shock  mounts. 

4.  Superimpose  a  reference  trace  on  e«ch  response  to  aid  in  the 
determination  of  the  base  line  and  in  indicating  any  drift  or  divergence 
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of  th«  r««poniie  tract. 

5.  InreBtlgate  hj   flight  t»et  the  lateral  ■otlon  perfonwnce 
functions  of  an  airplane  harln^  a  relatively  high  degree  of  epiral 
dlTergence.  Since  the  degree  of  spiral  diTergence  Increases  with  angle 
of  attack,  it  aay  be  necessary  to  perfom  these  measurements  at  a  high 
angle  of  attack  in  order  to  obtain  suitable  responses. 

6.  Inrestigate  the  performance  functions  of  an  airplane  haring 
dynaalcs  quite  different  frou  those  of  the  B-25J .  Sereral  uodom  filter 
types  exhibit  an  undaaped  oscillation  about  the  Z  axle.  A  study  of  this 
■otlon  should  prore  most  worthwhile  and  interesting. 
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APPElfDDC  A 
INSCRIPTION  07  AIBPLANI  IN  TIST  CONTIOUEATION 
The  airplane  used  In  the  tests  was  a  USAF  B-25J  under  bailment  contract 
to  the  In  at  rumen  tat  Ion  Laboratory,  M.I.T.  The  airplane  had  prerloualy  "been 
stripped  of  all  guns  and  turrets,  as  shown  In  ?lg.  A-1. 

Following  Is  a  listing  of  the  specif ications  and  dimensions  of  the  air- 
plane as  It  was  used  in  the  teste  described  in  this  thesis: 
Manufacturer 


Serial  Number 
Center-of-gravlty  range 

Forward  limit 

Aft  liAlt 

For  tests 
Gear  retraction  moment/lOOO 
Overall  length 
Height 
Wing,  airfoil  section 

Boot 

Tip 

Span 

Area  (Total) 

Area  (less  ailerons) 

Aspect  ratio 

Taper  ratio  (tip  chord/root  chord) 

Chord  at  root 


North  American  Ariation 
Corporation,  Inc. 

B-25J-25-N.C. 

44-30328 

16.5  percent  M.A.C. 
34.0  percent  M.A.C. 
27.4  percent  M.A.C. 
59.7    in- lbs 
53  ft   5.75  in 
16  ft  4.19  in 

NACA  23017 

NACA  4409  R 

67   ft   6.704   in 

609.8   sq  ft 

677.67   sq  ft 

7.48 

0.415 

12   ft   10.6   in 


Chord  near  tip    (380.75  from  fuselage)      5  ft   10.14   in 
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Projected  tip  chord 

Mean  aerodynAiulc  chord 

Length 

Distance  of  leading  edge  "back  of 

noee  reference  datura  line 

Incidence 

Hoot 

Tip 
Dihedral  (25  percent  line) 

Wing  center  section 

Ving  outer  section 

Sweep^ack  (leading  edge) 

Aileron 

Type 

Area 

Span 

Chord 

Inboard 
Outboard 

Trarel 

Up  (from  neutral) 
Down  (from  neutral) 

Tab  area 

Tab  span 

Tab  chord  (mean) 

Tab  travel 

Up  (from  aileron  trailing  edge) 
Down  (from  aileron  trailing  edge) 


Wing  flaps 
Type 


6  ft  4.257  in 

116.16  In 
212.68  in 


3  0*  30 


1 1 


0  21'  39 


4  38*  23" 
0°  12*  39 


•  * 


Area  (total) 


4°  12*   13*  • 


Internal  balance  unsealed 
32.13  s<i  ft 
137.25  In 


19.87  in 
12.45  in 


28" 
15« 

3.86  sq  ft 

57.911  in 

5.375  in 

12« 
12^ 


Tariable-slot 
trailing  edge 

75.8  sq  ft 
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Center  section  flaps 

Area  32.6  tq  ft 

Span  (nean)  127.5  in 

Travel  (down)  43 

Oafboard  flap 

Area  43.2  tq  ft 

Span  (mean)  214.5  in 

V  o 

Travel  (down;  46 

Horizontal  tail  surface 

Area  (including  elevators)  132.4  sq  ft 

Span  22  ft  2  in 

Mailniim  chord  7  ft  1.875  in 

Incidence  2 

Dihedral  '    None 

llevator  area  (two,  Including  tabs) 

Ilevator  span  (two) 

llevator  travel 

Up  (from  streanline 
with  stabilizer) 
Down  (from  streamline 
with  stabilizer) 

llevator  trim  tab  area  (total) 

llevator  trim  tab  span  (2) 

llevator  trim  tab  mean  chord 

llevator  trim  tab  travel 

Up  (from  elevator  trailing  edge) 
Down  (from  elevator  trailing  edge) 

Vertical  tail  surfaces 

Area  (total)  91  sq  ft 

Tin  area  (two)  47 .B  sq  ft 

Span  (mazimaa)  104.25  in 
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60.6 

sq  ft 

216.! 

>  in 

26" 

10^ 

4.22 

•q  ft 

57.2 

in 

7.87* 

>'in 

4* 

20° 

Badd«r  area  (two  including  tabs) 
Budder  span  (aaxlBua) 


Saddtr  trarel 

Right  (froB  etreaidine  vlth  fin) 
Left  (froB  ttreABillne  with  fin) 

Rodder  trim  tab  area  (total) 

Rudder  trim  tab  spaa 

Rudder  trlB  tab  mean  chord 


Rudder  trim  tab  traTel 

Ri^ht  (froB  rudder  trailing  edge)   12 
Left  (froa  rudder  trailing  edge)   12 


43.2 

•q 

ft 

101.1  is 

I 

20^ 
20« 

3.18 

■q 

ft 

36.4 

in 

6.3  ! 

Ln 

fuselage 

Maxlmm  width 

4  ft  8.5  In 

MaximuiD  height 

7  ft  4  in 

Length  (tip  of  nose 

to  tip  of 

tail) 

53  ft  5.75  in 

Ingine 

Vunber 

2 

Type 

Wright  Cyclone 

Designation 

a-2600-13  or  -29 

loBber  of  cylinders 

i 

14  -  double  row 

Oear  ratio 

16:9 

Supercharger  gear  ra 
Low 
High 

itio 

7.06tl 
10 .06:1 

Power  settings: 
Condition 

rpB 

2600 

2400 

Manifold  pressure  (In  of  H  ) 

Take-off 
NHP 

46                « 
40 

Cruise 
Approach 
Tor  tests 

2100 
2000 
1800 

29 
20 
24 

Propellers 

Manufacturer 

HaBllton  Standard 

Type 


Controllable  pitch 
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Iiiabtr  of  blad*! 
riaaater 
Blade  type 
Control  type 

Hu"b  type 

Pitch  eettln^ 
Low  (fine) 
High  (eoaree) 

Landing  gear  (estimated  data) 

Main  gear  assembly 

Weight  of  retractable  portion 
Radius  of  gyration  of  retractable 
portion  about  retraction  plrot 

Hose  gear  aeseably 

Weight  of  retractable  portion 
Radius  of  gyration  of  retractable 
portion  about  retraction  plTot 

Weight 

Moaents  of  Inertia  (eetlnated) 


hi 


3 

12  ft  7   In 

6359  A- 18 

Hydromatlc   full 
feathering 

23150-473 


22 
90^ 


600  lb   (each) 
48  in 

200  lb 
46  In 
25,637  lb 

63.000     slug  ft' 

4 

120,000     slug  ft' 

0 

-1930       slug  ft' 
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oo 
(\J 

O 


If) 
OJ 

I 

LU 

_J 
Q- 

cr 


CO 
LU 

I- 


4Z 


•810  704  (67-6  709  ) 


405.352 


4*38' 23"   DIHEDRAL  AT 
INTERSECTION  OF  CHORD 
PLANE   AND   LEADING 
EDGE  OF  WING  OF  CENTER 
SECTION 


r LEADING    EDGE  WING  OF 

\       OUTER     PANEL    AT  WING   JOINT 


0*2r'J9"NEG   DIHEDRAL 
AT   INTERSECTION  OF 
CHORD  PLANE   8  LEADING 
EDGE    OF  OUTER  PANEL 


LEADING  EDGE  OF 
CENTER  SECTION 
WING  AT  4.  AIRPLANE 


-WING  REF  LINE 
AT  t    SHIP 


232" 


650" 


-*! 


648-l3/l6"(54-l}/i6') 


LEW  AT  t  SHIP-*- 
-  224' 


-FRONT  't    SPAR 


483a-- 


--4I84' 


-145-1/2"- 


t 


74-1/4 


k*|37f 


-NORMAL   GROSS  WEIGHT  C  G  ^<,6-5/B  U- 

LANDINC   GEAR  UP  !  T 

20 

-EXTREME   REARWARD  POSITION 
OF  C  G    LANDING   GEAR  UP 


104-1/4' 


T 


169,- 
«^-'^'^"   (l5'-9,{ 


30"  SMOOTH  CONTOUR^ 
3" 


STA  0 


-195 


STATIC    GROUND   LINE 


-NORMAL   GROSS  WEIGHT  C  G  -^ '^       j" 

LANDING   GEAR    DOWN  I 


EXTREME    FORWARD   POSITION 
OF  C  G     LANDING   GEAR    DOWN 

47'   SMOOTH    CONTOUR 


-t   RUDDER    HINGE 


-13*44 


FIG.  A-2.     TEST    VEHICLE    DIMENSIONS 
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APFXIDIX  B 
IISTHUMETTATIOH  AFD  CALIBRATIOH 
A.   IHSTHUWCNTATiai 

As  aoted  in  the  Introduction,  finding  the  pcrfornanee  function  by  aeant 
of  a  pulse  input  has  the  adranta^es  of  keeping  the  response  within  the 
linearity  linlts  of  deriation  from  the  equilihriuji  flight  conditions,  and  of 
requiring  a  minioua  of  instruaentation.  The  derices  required  for  applying 
the  input  pulses  were  designed  to  accomplish  the  following  functions: 

1.  Allow  rotation  of  the  control  surfaces  between  any  desired  limits* 
with  positire  stops  at  the  ends  of  the  arc  of  trarel. 

2.  Allow  rariations  in  the  size  of  the  pulse. 

3.  Allow  the  pulse  to  he  applied  in  both  directions}  that  is,  allow 
changes  in  the  position  of  the  mechanism  relatire  to  the  pilot's 
controls  without  altering  the  pulse  size.  This  function  also 
allows  the  pilot  to  hold  two  controls  rigidly  against  stops  while 
the  aircraft  is  responding  to  a  pulse  imparted  by  the  third  control. 

4.  Allow  quick  disconnect  so  as  to  make  it  possible  to  remore  all 
impediments  from  full  use  of  the  pilot's  controls  in  case  of 
emergency. 

The  mechanisms  employed  were  chosen  from  many  possible  designs  as  the 

simplest  that  could  perform  the  above  functions  and  still  be  fabricated  with 

a  minimum  of  machining.  The  aileron  pulsing  mechanism,  shown  in  Figures  B-1, 

B-2,  and  B-3,  consists  of  a  plate  mounted  on  the  right  control  wheel  which  is 

slotted  to  receire  two  piroted  bars  which  are  set  to  gire  the  desired  sise  of 

pulse.  A  slotted  length  of  duralumin  angle,  hinged  at  one  edge,  was  mounted 

on  a  shelf  fastened  to  the  top  of  thA  control  colxuon.  The  slot  allowed  a 

single  stop  to  be  positioned  so  that  the  chosen  pulse  could  be  initiated  at 

any  trim  condition  by  moving  the  single  stop  against  either  one  of  the  stops 

on  the  wheel  plate.  The  hinged  edge  allows  the  stop  to  be  thrown  out  of 

position  in  order  to  clear  the  aileron  control  in  case  of  emergency  in  flight. 

Approximately  45  degrees  of  wheel  rotation  were  used  in  applying  the  desired 

aileron  pulse* 


Th«  ■•chAnltB  used  for  applying  rudder  pulees  It  shown  in  Tlguree 
B-4  and  B-5,  Choice  of  the  design  selected  was  dictated  by  the  limited 
installation  space  and  by  the  requirement  that  the  pulse  setting  and  quick 
disconnect  units  must  "be  easily  accessible  to  the  co-pilot  in  flight.  The 
laechanlsin  in  its  final  configuration  consists  of  a  long  rod,  one  end  of  which 
is  attached  to  the  right  rudder  pedal.  The  other  end  of  the  rod  slides  in  a 
block  which  is  plToted  at  the  center,  and  mounted  on  a  plate  fixed  to  the 
aircraft  structure.  The  block  receirlng  the  sliding  rod  is  cut  away  at  the 
top  80  that  a  stop  pin  welded  to  the  sliding  rod  may  move  in  the  slot  thus 
provided .   The  cover  of  the  block  is  made  up  of  a  frame  hinged  at  the  inboard 
edge  and  containing  two  stop  bars  connected  by  a  positioning  screw  to  provide 
adjustment  of  pulse  size.   In  addition,  both  bars,  holding  the  set  pulse  size, 
may  be  moved  by  means  of  a  second  screw  to  permit  engaging  the  pin  on  the 
sliding  rod  with  either  of  the  stop  bars  when  in  the  trim  condition.  The 
hinge  on  the  inboard  side  allows  the  block  cover  to  be  moved  clear  of  the  rod 
stop,  thus  quickly  and  positively  freeing  the  rudder  controls.  Figures  B-4 
and  B-6  show  the  cover  block  clear  of  the  rod  stop  In  the  stowed  condition. 
Approximately  Ci3  in.  of  linear  throw  of  the  rudder  pedals  provided  the 
desired  rudder  pulse  size. 

The  mlcrosyn  pickoff  units  installed  for  measurement  of  control  surface 
rotation  were  of  the  variable  transformer  type  designed  by  Dr.  fi.I.  Mueller 
of  the  M.I.T.  Instrumentation  Laboratory.   The  Installation  of  these  xinlti 
in  the  right  nnd  left  wing  are  shown  in  Figures  B-6  and  B-7  respectively. 
The  mounting  positions  and  linkage  lengths  are  Identical  in  each  wing.  The 
final  determination  of  linkage  sensitivity,  that  is,  piokoff  rotation  per 
control  surface  rotation,  was  governed  by  accessibility  in  the  wing  or  rudder 
cavity,  linearity  limits  of  the  pickoff  (approximately  150  mlllirradians)  ,«Lnd  the 
maximum  control  surface  motion  anticipated  in  obtaining  the  desired  magnitude 
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of  reeponsc.   Shielded  plckoff  leads  were  used  throughout,  and,  p.e  ehovn 
In  Tig9  .   B-6  and  B-7 ,  safety  braclcett  were  attached  to  the  mlcroeyTi  BOuntlngB 
to  prevent  any  possible  damage  or  Jamming  of  controls  in  case  th*>  linkage 
worked  free  at  the  pickoff  rotor. 

The  rudder  plckoff  Installation  is  shown  in  Tig.  B-B .   The  linkage  is 

t 

similar  to  those  used  in  the  win^s,  the  bell  crank  at  the  radder  control 

serving  as  the  initial  point  of  motion  measurement. 

The  rate  gyros  employed  were  elevation  and  deflection  gyros  developed 
by  the  Instrumentation  Laboratory,  M.I.T.   The  characteristics  of  these 
gyros  were  modified  to  suit  the  measurement  requirements  of  the  project  by 
changing  the  damping,  and  were  installed  as  roll  and  yaw  rate   gyros 
respectively. 

Laboratory  calibration  of  the  gyros  prior  to  installfttion  In  the  aircraft 
was  necessary  to  determine  the  excitation  currents  desired  for  the  elastic 
restraints  and  the  signal  generators  (pick-off s),  the  tenperatupt  required 
for  the  damping  fluid,  and  the  friction  torque  level. 

The  friction   torque  level  was  first  determined  to  Judge  whether  the  gyros 
would  provide  a  satisfactorily  sensitive  response  to  allow  measurements  of 
the  desired  accuracy  and  precision.   The  final  results  are  shown  in  Tiguree  9 
and  10  for  the  roll  and  yaw  rate  gyros,  respectively.   The  maximum  friction 
torque  level  for  the  roll  gyro  was  50  dy  cm,  and  for  the  yaw  gyro  was  140 
dy  cms,  and  the  unbalance  was  reduced  to  a  torque  level  below  that  of  the 
friction  level.  The  maximum  friction  torque  levrl  for  th#  roll  rate  gyro 
represents  an  angular  velocity  input  of  0.25   mllli-radians/soc  and  for  the  yaw 
rate  gyro,  the  maximum  friction  torque  represents  an  angular  input  velocity  of 
0.7  milli-radians/sec .  This  friction  level  is  considered  satisfactory  since 
the  oscillograph  record  ordinates  can  be  measured  only  to  a  maxlmuir.  accuracy  of 
.01  inch,  the  latter  figure  representing  approximately  1/2  milli-radlans/sec  in 
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■oet  of  th«  recordt. 

The  rate  gyro  equations  of  motion  were  Investigated  to  malce  certain 
that  the  desired  characteristic  could  he  obtained.  The  general  equation 
representing  both  gyros  is: 

•  •  • 

*(  )  ^^llnk)(M  M  )  "  ^(c)(eff)  ^  "  °c  ^c  "  ^6r)(A  H  )  \ 
g  c  c  c 

where  the  subscript  c  refers  the  indicated  quantity  to  the  computer  shaft, 
and  the  subscript  g  refers  the  indicated  quantity  to  the  glmbal  shaft. 

This  equation  may  be  reduced   to  the  general  form  of  the  second-order 
differential  equation  yielding 

^  s,  -^ i 


e 
e 

^^c(eff)  ^n 

V          ^^- 

^n  •        ^ 

^er)(AM) 
^c(.ff) 

where 

(DH)    : 

(FE)   . 

To  determine  the  gyro  undairped  natural  frequency,  the  effective  moment 
of  inertia  at  the  computer  shaft  was  measured  using  low  elastic  restraint 
currents  aiul  no  danping.  The  period  was  recorded  on  a  brush  type  oscillograph, 
and  the  monent  of  inertia  was  then  determined  from  the  following  equation: 

T         ^(er)(AM)  -2    ,-.  .         2 
^c(eff)  "      '2 ^  -  264  dyne-cirsec 

Hechecke  at  different  values  of  elastic  restraint  verified  the  above 
figure.  To  establish  the  amount  of  elastic  restraint  desired  on  the  gyros, 
mailmum  rates  of  roll  and  yaw  for  the  magnitude  of  control  surface  deflection 


_eb. 

2U 
Cornell  Laboratory  Reports  covering  fate  responses  to  step  inputs  of  control 


inputs  anticipated  were  determined  from  B-26  roll  rate  data  (  L,  )  and  from 
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deflection  (Befe.  4  and  5).  The  maxinan  rates  deternlned  vere  0.4  rad/eec 
in  roll  and  0.5  rad/seo  in  jav.  The  elastic  restraint  unite  were  then 
calitrated  for  "both  gyros,  ae  shown  in  Tlgg.  B-11  and  B-12,  using  various 
ralues  of  elastic  restraint  cvirrent.  The  maziruB  rate  tabulated  for  each 
sensitiTity  indicates  the  rate  input  for  the  particular  sensitlrity  which 
will  rotate  the  gyro  computer  shaft  150  milliradiant.  At  deflections 
exceeding  this  value,  the  signal  generators  "become  nonlinear.  The  ovei-alj 
sensitivities  of  the  elastic  restraints,  considering  both  excitation  current 
and  angulf^r  displacement  as  inputs  are: 

(rol])(er)(A,  1  ;  M)        (^j2  ^ 

Selection  of  the  desired  undamped  natural  frequency  depends  on  the 
desired  frequency  response  of  the  gyro.  The  rate  gyros  aast  accurately 
reproduce  aircraft  response  below  about  10  rads/sec,  offering  negligible 
or  BinlBUK  attenuation  and  phase  shift*  while  at  the  resonant  structural 
vibration  frequencies  of  the  airplane,  the  attenuation  Bust  be  a  aaziiaaB. 
Previous  instrumentation  work  on  the  B-25  has  Indicated  that  the  most 
pronounced  resonant  frequencies  were  fooind  at  16  and  30  cps.  To  fulfill 
the  response  requirements  outlined  above,  the  following  acceptable  limits 
for  the  gyro  characterietics  were  determiced: 

Undamped  Natural  Trequency,  n     ■  B  -  8  cps 
Damping  Ratio,  (DH)  -  0.6  -  1.0 
With  the  above  ranges,  elastic  restraint  ourrente  were  selected  to  satisfy 
both  undamped  natural  frequency  emd  maxinuB  rate  input  requirements. 

A  damping  fluid  was  then  selected  to  fit  the  abeve  requirements  of  DS. 
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Tigur*  B-13  thowB  the  rlscoBlty  tb  temperature  chamcterlBticB  of  damping 
fluid  DCC  350,  which  wsb  the  fluid  uaed  In  "both  gyroB .  DCC  350  offered  the 
Bmalleet  gmdient  of  TlBCoeitj  to  temperature  of  any  of  the  fluids  available 
haring  the  rlBCOsity  range  desired. 

The  following  chracterlBtica ,  coefficients,  and  sensitlTltles  describt 
the  gyros  as  used  in  the  final  test  flights: 

Roll  rate  gyro       Yaw  rate  fnyro 

2  2 

1/  w  --\  264  dy  cm  sec        264  dy  cm  sec 

Of   ^  17200 -Jf^^       24000^/°' 

(c;  rad/sec  rad/sec 

0  _  0 

Daaper  tenperatur*  160  J  177  I 

8/  ^f.u^  3.8X10^^^      4.66X10^^^ 

(•r)(AM)  r&d  rad 

V  38  rad/seo  42  rad/sec 

n 

B  6,06  cps  6,7  cps 

n 

(DE)  .857  l.Oe 

The  gyros  were  mounted  in  a  dustproof  container  and  secured  in  the  air- 
plane on  an  adjustable  mount  which  allowed  aligning  the  gyros  with  the 
horizontal.  Since  the  B-25  trims  at  175  mph  in  a  nose  high  attitude  with 
its  horizontal  reference  line  3.6  above  the  horizontal,  the  gyro  container 
must  "be  tilted  accordingly.  The  adjustable  mount  is  shown  in  7lg.  B-14, 
and  the  gyro  container  in  position  on  the  mount  is  shown  in  Tig.  B-15.  The 
gyro  container  utilized  three  thermostat ically  controlled  heaters  to  maintain 
a  fixed  a»blent  temperature  for  the  gyros, 

Tigure  B-16  shows  the  amplifier  panel.  Four  channels  of  amplification 
were  used  to  amplify  the  pict-off  voltages  representing  total  aileron  deflection, 
rudder  deflection,  angular  velocity  in  yaw,  and  angular  velocity  in  roll, 
lach  amplifier  channel  performed  three  functions  (l)  amplified  the  AC  picfc-off 
signal,  (2)  demodulated  the  AC  signal,  and  (3)  amplified  the  resulting  DC 

signal.  The  output  of  each  amplifier  channel  la  controlled  by  12  attenuation 
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■•tiln^t  In  order  to  adjust  oscillograph  current  to  gire  the  desired 
oscillograph  displacenent . 

The  oscillograph  used  is  shown  in  ^ig.  B-17 .   It  is  a  12-channel 
Consolidated  Recording  Oscillograph.   The  dynamics  of  the  oscillograph 
galTemoueter  units  hare  a  nefligi'ble  distortion  effect  In  reproducing  the 
Inputs  introduced  tj   the  amplifier  outputs.  The  natural  frequency  of  the 
units  used  is  120  cpt  and  the  DB  is  0.7  approz. 

The  stiffness  current  regulator  shovn  mounted  in  position  in  Tig.lrld 
Is  a  closed-loop  feedback  regulator  functioning  to  hold  the  stiffness  current 
at  a  set  current  within  1  percent  Tariation.  figures  B-IB  and  B-17  show  the 
exposed  regulator  chassis,  the  chassis  haring  been  remoTed  from  its  case. 

To  integrate  the  added  instrumentation  with  the  instrumentation 
prerioualy  existing  in  the  aircraft,  an  additional  Junction  box  was 
Installed.   It's  location  is  noted  in  Tig.  B-17. 

The  three-phase,  400-cycle  thyratron  controlled  inverter  initially 
installed  was  found  to  give  unsatisfactory  regulation  for  the  needs  of 
this  project.  Consequently ,  It  was  replaced  by  a  USAT  Ptf-16  single-phase 
inyerter,  which  was  found  to  give  satisfactory  regulation.  This  change 
necessitated  use  of  a  phase  splitter  and  stepdown  transformer  to  obtain 
the  3  phase,  28y  gyro  wheel  supply.  A  block  diagrams  of  the  overall 
instrumentation  system  is  shown  in  Figure  B-26 . 

As  shown  In  Fig.  B-26  the  output  of  each  piokoff  was  loaded  with  a 
1000-ohm  resistance  to  reduce  the  noise  level  and  the  null  voltage.  A 
variable  resistance  and  a  variable  inductance  were  placed  in  series  with 
the  pickoff  excitation  coils  (P.O.  coils  in  series).  The  resistance  was 
used  to  set  the  desired  excitation  current,  and  the  inductance  was  used  to 
bring  the  pickoff  outputs  into  proper  phase. 
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B.   CALIBRATION 

The  objective  of  the  mcRBurlng  Bysteni  calibration  for  thi«  project 
vae  to  eetablish  orerall  BeneltlvltieB,  expressing  Inputs  to  the  measuring 
Instnments  In  terms  of  output  linear  displacements  at  the  osclllo<»raphi . 
The  senBltlvities  are  expressed  in  such  a  manner  that  multiplication  of  any 
displacement  on  the  osclllogmph  record  by  the  appropriate  sensitirity  yields 
the  measured  quantity  In  the  desired  uniti.  These  eensltlTltles  are: 

^(all)(X   a  )         '^•^^^^ 

osc  A 
Sr^^Wr    A  ^  deg/ln 


(rttd)(X    C  ) 
osc  r 


S 


deg/sec 


(roll  gyTo)(X_,  W  )      in 

g  deg/see 

(yaw  gyTo)(l^_^  W  )        In 

osc  it 

For  purposes  of  Bimplicity  and  system  flexibility,  the  overall  sensitivities 
listed  above  were  broken  Into  two  component  sensitivities,  one  representing 
the  Instrument  sensitivity  and  one  representing  the  Indicating  system 
sensitivity.  The  Instrument  sensitivity  is  a  fixed  quantity  determined  In 
the  case  of  the  control  surface  rotations  by  the  llnka^jje  sensitivities  and 
the  signal  generator  •xcitatlon  current,  and  in  the  case  of  the  rate  gyros 
by  the  elastic  restraint  current  and  the  signal  generator  •■citation  current. 

The  Indicating  system  sensitivity  may  be  varied  at  will  by  changing  the 
attenuation  setting  at  the  amplifier.   The  amplifier  sensitivity  allows  the 
desired  size  of  recorded  input  pulses  and  output  angular  velocities  to  be 
selected.   Selection  of  the  proper  size  of  recorded  dlsplpcements  is  governed 
on  one  hand  by  the  linearity  limits  of  the  oscillof^raph  trace  and  on  the 
other  hand  by  the  size  record  required  for  ease  of  analysis. 

In  order  to  obviate  the  necessity  for  calibration  of  each  amplifier  at 
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•Ach  attenuation   tettln^,   a    * '■fttter* *   pickoff  and  a  calibration  Tolta«« 

•oui*c«,    the   latter  prorldln^  12  roltftge   stopa,   wore   employed.      The   naet«r 

pickoff  was  calibrated  (Ti^re  B-18)  to  deterulne  the  llmlte  of  linearity. 

The  •ensltlTlty  determined  was  S,   w.        x  «  4  T,^?^?*   ,  with  300 

(■poMA    •   ;     division 
npo  inpo 

Baeter  -p.o.  dlrltlons  as  the  limit  of  the  linearity  ran^e.  The  master 
pickoff  was  then  used  to  * '"buck-out  *  *  or  neutralize  the  voltage  produced 
"by  the  instrument  plokefft.  Thus*  instrument  sensltlrlties  vera  obtained 
both  in  terms  of  ralta^e  output  for  Input  rotation  or  ani^ular  relocity  and 
in  terms  of  master  pickoff  output  an^le  for  input  rotation  or  angular  relocity. 

The  master  pickoff  Is  shown  in  Fig.  B-19 .  The  master  pickoff  could  be 
plugged  into  aoiy  of  the  pickoff  outputs  by  means  of  bayonet  Jack  receptlcle^ 
on  the  amplifier  panel.  As  shown  in  Fig.  B-26,  the  particular  pickoff  oat- 
put  could  then  be  switched  from  the  amplifiers  to  the  master  pickoff. 

The  procedure  for  obtaining  the  instrument  sensltlTities  vas  ftraii^t- 
forward.  For  the  ailerons,  a  propeller  bubble  protractor,  measuring  to  O.l 
degree  accuracy  was  used  to  measure  the  angle  of  both  ailerons  with  respect 
to  the  horizontal,  and  the  total  aileron  differential  angle  thus  obtained. 
For  each  measured  6  ,  the  corresponding  roltage  was  recorded.  This  roltage 
was  then  nulled  by  the  MPO  and  the  corresponding  angle  of  the  MPO  recorded. 
The  Toltage  representing  total  aileron  angle  was  achlered  by  adding  the 
Yoltages  of  the  two  aileron  pickoff  outputs  in  series. 

For  calibration  of  the  rudder  deflection  to  MPO  angle,  the  same  system 
iras  used  as  outlined  for  the  ailerons  with  the  exception  of  the  means  of 
measuring  rudder  deflection.  Figure  B-20  shows  the  method  employed.  The 
center  of  rotation  of  the  rudder  was  determined  by  means  of  a  plumb  bob. 
A  piece  of  mirror  was  then  attached  to  the  rudder  at  the  center  of  rotation, 
^nd  the  plumb  bob  lined  up  with  the  center  of  a  sector  of  a  circle,  the 
latter  being  graduated  in  mils.  Measurements  were  then  read  by  means  of  a 


telescope  equipped  with  a  hair  line, 

rigure  B-21  shows  the  set-up  used  to  calibrate  the  gjros .   On  a  small 

turntable,  having  a  selection  of  4  speeds,  an  Inclined  plate  was  mounted  to 

receive  the  gyro  container.   The  can  was  so  aligned  that  the  spin  axis  of 

the  roll  rate  gyro  was  rotated  a  measured  angle  about  the  gimbal  ails.   With 

the  two  gyros  properly  aligned  within  the  can,  this  also  rotated  the  input 

ajcis  of  the  yaw  gyro  about  its  spin  axis  by  the  same  measured  angle.  As  a 

consequenct,  the  input  to  each  gyro  was  calculated  in  the  following  manner: 

Holl  gyro       Wj  -  W^^  sin  6 

Taw  gyro        1^  =  *jj  cos  « 

where  W   if  the  angular  Telocity 

of  the  turntable,  and  6  is  the  measured  angle  of  tilt. 

The  calibration  curves  for  all  four  measuring  Instruuents  showing  the 
measured  quantity  vs  MPO  divisions  (each  div.  =  approx.  2  milt)  and,  for 
the  gyros,  the  measured  quantity  vs  plcfcoff  output  voltage  are  shown  in 
Tigs.  3-22  through  B-25. 

The  Instrument  senaitlvitiea  thus  determined  are 

[aileron]   ^^^^l^,^^^  tj   ^(MPO)(a^j^  »^)  =  ^all)(M)(i^  S^)  "  °-°2"  dlf 
[nxdd.r]    S(^)(^^  ,^)  ^MPO)(A^  .^)  '   ^rud)(M)(A^  «^)  "  °-°°«^«  flf 
[roll  gyro]  S(^^j(^^  ^^j  '(HPO)(A^  .^)  •  ^(rg)(M)(A^  Wj)  ■  °-5«^  Til'i^" 

[7a.  gyro]  Sj^^j^^^  ^^)  '(«PO)(A^  .„)  =  '(rg)(V)(A^  Wj)   '•'''  ^^  «P0 

The  indicating  system  sensitivities  were  determined  on  each  flight 
whenever  data  was  recorded.   The  indicating  eystera  sensitivity  is  made  up 
of  the  amplifier  sensitlTlty,  oscillograph  sensitivity  and  master  plcfcoff 
sensitivity,  that  is: 
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Upon  the  completion  of  each  eerleB  of  flight  rune  the  meaaurln^ 
Inatnifflents  vere  switched  Rway  from  the  inputs  to  the  aarpllflere,  and  a 
12-8tep  voltage  eoiirce  switched  onto  the  aaplifier  Inpute   (shown  in  Ti^.B-19)  , 
As  the  Toltfl^e  was  rarled  through  the  12  steps,  a  record  was  made  on  the 
oscillograph.  These  same  12  rolta^e  steps  were  then  applied  to  the  master 
pickoff,  the  latter  instrument  beln^  rotated  until  each  calibration  rolta^e 
was  ntllled  and  the  corresponding  output  angle  of  the  master  piclcoff  was 
recorded.  Thus  increments  of  oscillograph  displacement  corresponding  to  in- 
crementt  of  MPO  angle  were  available  for  each  test  fli^t ,  and  an  OTersdl 
Indicating  system  sensitirity  was  arailable  without  knowledge  of  the 
oscillograph  and  sunplifier  component  sensitirities .   That  is, 

The  orerall  system  sensitirity  for  each  measuring  channel  was  thus 
arailable  by  eonbiBin^  instruaent  and  indicating  system  sensitiritiet . 
Therefore,  the  general  expression  of  the  sensitirity  is 

»(,  eh«.)(l^,^  ,)  ■  »(  ){M)(1^,  ,)  \   )(lnd  .ye)(X„^  i^) 

To  determine  the  system  errors  and  anoertaintiet,  static  tests  were 
made  on  the  ground.  A  known  ralue  of  each  quantity  to  be  measured  was 
applied  to  the  respective  measuring  instruments,  and  an  oscillograph 
record  made  of  the  outputs.  The  calibration  source  voltages  were  then 
recorded  at  the  same  attenuation  settings,  and  a  record  of  master  pickoff 
angles  recorded  for  each  calibration  voltage.   In  this  manner  the  overall 
sensitivity  could  be  compared  with  the  product  of  the  instrument  sensitivity 
and  the  indicating  system  sensitivity. 
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The  linearity  limite  arvd  of  the  meaeurln^  system  components  are  ab 
follows: 

Component  Linearity  Hange 

MMbsr  P.O.  ♦SOO  dirs  (160  ails) 
Aileron  Defl,       30  total  deflection 

Hudder  Defl.  ^7  rudder  deflection 

Roll  gyro  ^0.3  rad/sec 

Taw  gyro  ^.35  rad/sec 

Oscillograph  defl.  ^S.?  in  (approx.) 

The  largest  source  of  error  in  the  measuring  system  lies  in  the 
inseneitirity  of  the  Yoltmeter  used  as  an  indicator  when  nulling  a  roltage 
with  the  master  p.o.  If forts  to  reproduce  results  have  produced  rariations 
as  high  as  4  dirisions  of  the  MPO,  with  an  average  error  of  1  to  1  l/2 
divisions.  This  represents  am  enror,  or  null  voltage,  of  4  to  6  bt  in  matching 
voltages. 

A  change  in  the  instrumentation  was  effected  Just  prior  to  the  final 
test  flight  (See  Flight  Record  -  Appendix  D) .  The  gyro  container  was 
secured  rigidly  to  the  aircraft  structure  during  flight  in  order  to  eliminate 
recording  vibrations  ittip  "by  the  gyro  container  mount  and  the  Lord  shock 
mounts.  This  modification  resulted  in  removing  an  undesirable  10  cps 
transient  from  the  records  which  occurred  whenever  the  aircraft  was  respond- 
ing to  input  pulses.  This  change  made  it  possible  to  asJoe  the  shock  mounts 
ineffective  during  flight,  and  yet  to  free  them  when  desired  in  order  to 
protect  the  gyros  from  landing  shock. 
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F,G   BI4    ADJUSTABLE    MOUNT  FOR    GYRO   C^^NTAINER 
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IFPXMDIX  C 
ANALYSIS  07  A  DITIsaiNT  TRANS IINT  EXSPONSX 
Reference*  (2)  and  (3)  point  out  that  the  approximation  of  the  Toorler 
Tranefonn  is  applicable  only  to  stable  eystent.  This  car  be  easily  seen  frou 
a  consideration  of  the  transient  response  of  a  syetem  described  by  a 
characteristic  equation  harin^  a  poeitlTS  root.  Tor  such  a  reeponee,  the  sux 

nBi,<;,.,, 

will  not  converge.  This  condition  is  found  in  the  lateral  notion  of  aany 
airplanes  I  where  the  positlTe  root  is  aspoclated  with  the  motion  called 
spiral  dlverjtence. 

Since  the  B-26J  used  in  the  present  inrestigatlon  did  not  exhibit  thie 
type  of  response,  a  simple  dynamic  system  illustrating  divergent  motion  is 
analyzed  In  this  appendix,  using  the  finite  pulse  method. 

Consider  the  system  seheBatically  pictured  in  Tig.  C-l.  The  Base,  n, 
in  this  figure  is  assumed  to  be  supported  in  the  gravity  field  by  a  rigid, 
messless, rod.  The  Laplace  Transform  of  the  characteristic  equation  for  the 
^  small  angles  of  motion  considered  in  Fig.  C-1,  (Ref.  8),  is 
[b^8^  -^  C>^  8  -  Big]  0  =0  (C-1) 

Assume  that  a  forcing  function  is  applied  consisting  of  a  unit  triangle 
of  force,  UT(t  -  At)  which  is  described  by 

T  -  0  0  >  t 

T  »  -^  t  0  4  t  ^  At 

At 

I*2-rA-t  AT^ti2AT 

At 

T  -  0  t  >  2At 

At  "  0.1  second 

-83- 


The  transformed  equation  of  motion  "becomeB 

[mi  t^  -f  C  i  •  -  iB«]»(«Jfc=  (LT)[UT(t  -  At)]     (C-2) 

where 

(LT)[UT(t  -  At)]=  ^ 5-^^ 

At  n 


Solving  for  0  (§) , 
®  (•)  = 


(1  -  g    ) 


(At   ,^)(ni)(.^  ^^  b  -■!-) 


(0-3) 


In  thle  inveetigation  of  the  B-25J  lateral  motion,  the  reeponee  which 

was  measured  was  the  angular  rate  In  roll  or  yaw.  Tor  this  reason,  the 

time  rate  of  change  of  the  angular  position®  will  be  considered  as  the  output 

in  thle  example.  The  rate  of  change  of  this  output,  ®  ,  is  glren  by 

/,    -Ats\2 

®  (s)  =    '^^"  ^  ^  ^ (C-4) 

At  tti  tin     *—•  -;5-) 

factoring  the  characteristic  equation, 


Let 


5-  J  =1  *  i- 


2m     „  ^ 


b*t 


Bote  that  a  is  always  poaltire  and    t^iat    a  and  b  are  poeltlTe,  real 

numbers* 

By  appropriate  substitution  of  a  and  b,   the  eqxiatlon  of  motion,   eq .    (C-4) 

can  be  written, 

/ ,  -At  8  N  2 

0(e)    =     ^^-^^-5 ^ (0-5) 

At     mi    s    (s  -   a)(B  •*■  b) 
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The  perfornaac©  function  (Pr)r_  A  -i  1«  d«fln*d  ae 


•  -  J» 


Then, 


A  ralue  of  0.1  vat  aeeigned  to  a,  and  a  ralue  of  2.0  to  b.   The  ralae 
of  0.1  for  a  1«  such  that  the  antplltude  of  the  dirergence  will  doable  In  6.93 
0econdt.  Some  airplanes  hare  this  ouch  dlTorgence  at  hl^  anklet  of  attack. 
The  Talne  of  2.0  for  h  Is  somewhat  snail  compared  with  that  for  the  highly 
damped  motion  of  most  airplanes.  The  major  consideration  in  the  eholee  of 
these  roots  was  to  choose  valaes  that  wonld  reduce  the  calculations  as  waeh 
ae  praetieahle,  and  yet  that  would  retain  some  ••■blanoe  to  those  encountered 
in  airplane  lateral  performance  functions* 

These  choices  were  consistent  with  the  choice  of  n  and  /  as  1  and  2 
respect irely.  The  units  are  those  of  any  consistent  set  of  physical  units. 

Substituting  these  numbers  and  letting 
e  «  J» 
the  performance  function  becomes 

^^\f,   0  ]  "  F  (JCD  -  O.DO  +  2)         (C-7) 


This  performance  function  is  plotted  in  Tig.  C-6. 

Returning  to  eq .  (C-5) ,  by  the  inrerse  Laplace  transformation  the  time 
response  was  obtained  and  plotted  in  Tig.  C-2. 

In  Tig.  C-3,  the  portion  of  the  curre  of  Tig.  C-2  beyond  B  seconds  was 
plotted  on  semi-log  paper,  thus  enabling  determination  of  the  dlrergent  root 
The  straight  line  was  extended  to  zero  time  and  the  Intercept  found  .  Haying 
thus  found  an  expression  for  the  dlyergent  component  of  the  response,  It  was 
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•ubtracted  fron  tho  total  reeponte  learlng  the  remainder  shovn  in  7lg.  C-4. 
The  Integral  of  thin  remainder  is  conrergent  and  its  approximate  Tourier 
Tranefonu  ii 

^^•"^  ^   n  .  1.2...  ° 

fince  the  system  is  linear,  the  performance  function  determined  fro« 
the  response  cunre  can  "be  expressed  as  the  sum  of  components,  i.e., 

The  last  two  terms  on  the  rifht  side  of  eq.  (C-8)  were  determined  bj 
approximating  the  Fourier  Transform  "by  triangles.  The  first  term  can  be 
determined  since  the  analytical  expression  for  4.^(^)^4      »  !■  Icnovn,  and 
sine* 

aT)[q^^(t)]  =  (rT)[q^^(t)] 

This  analytical  expression  for  the  dirergent  component  of  the  response  as 
determined  from  the  time  response  is 

and 


I 


C».  9](dlT)    (j„-o.l)(»T)[OT(t)].-J^ 
From  this,  the  amplitude  ratio  eoid  phase  angle  can  be  calculated  for  aay 
given  •  . 

The  desired  performance 'function  can  be  obtained  by  performing  the 
addition  of  complex  quantities  indicated  in  eq.  (C-8) , 

Figure  C-5  shows  the  comparison  of  the  original  perfonnance  function 

with  the  performance  function  obtained  by  analysis  of  the  time  response.  The 
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Acouraoy  with  whioh  th«  original  p«rforMuic«  function  was  daplicAted  shows 
that  this  is  a  us«a^lt  method  of  deterainin^  a  p«rfomance  functibn  fron  a 
dirergent   time  response* 
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///////////////////////////////////////^^^ 


FIG.   C-1    SCHEMATIC   DIAGRAM    OF    SIMPLE    DYNAMICAL    SYSTEM 
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AFFUmiX  I) 


DATA  HIDUCTION  AHD  ITALUATICW 

The  method  used  to  obtain  a  deelred  perfomance  function,  knowing  the 
tlm«  response  of  a  aystta  to  a  known  Input  pulse,  is  outlined  In  Chapter  I. 
In  that  chapter  the  Method  of  approxlraatlnij  the  aathenatlcal  expressions 
graphlcallj  (within  the  desired  degree  of  accuracy)  Is  coTered.  Reference  3 
shows  the  method  applied  to  a  critically  damped  second  order  systen.  and  out- 
lines a  method  of  attack  to  "be  used  when  more  lightly  damped  oscillations 
are  present  in  the  time  response.  The  process  as  abore  outlined  has  been 
applied  to  the  oscillograph  test  records  obtained  In  this  inrestlgation. 

The  oscillograph  records  prorlded  the  following  Information: 


Input  Puleee 
Aileron  deflection,  «^  (♦  defl) 
Aileron  deflection,  6  (-  defl) 
Rudder  deflection,  «  (•♦■  defl) 

Rudder  deflection,  8^  (-  defl) 

r 


Output  Rest)onie 
Angular  rel.  In  yaw  and  roll 
Angular  rel.  in  yaw  and  roll 
Angular  rel.  In  yaw  and  roll 
Angular  Tel.  in  yaw  and  roll 


Sign  conrentions  used  in  designating  pulses  correspond  to  HACA  usagej  i.e., 
(•►)  Aileron  defl.  -  right  aileron  depressed 
(♦)  Rudder  defl.  -  rudder  trailii-g  ed^e  moved  to  the  left. 
Redaction  of  the  abore  listed  recorded  data  would  yield  the  following 
aircraft  performance  functions: 

(1)  Cw]^)(,  ^,  Wj) 


<2>  t"](A)(«  .  »,) 
a  t* 


(3)  W(^)(j^^,^) 


(«)  tW](i)(.^.  Wj) 
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The  tisM  crailable  for  redaction  of  data  was  not  sufficient  to  allow 


all  these  performance  functions  to  he  determined.  Therefore,  since  this 
project  is  prinarily  concerned  with  the  application  of  the  pulse  technique 
AS  a  means  of  determining  the  performance  function  of  any  dynamical  system, 
and  since  knowledge  of  the  actual  dynamics  of  the  test  rehicle  is  of 
secondary  importance,  the  records  were  analysed  to  achiere  the  following 
ohjectirest 

(1)  To  determine  the  degree  of  accuracy  or  uncertainty  in  repeating 
results. 

(2)  To  compare  the  results  of  determining  the  performance  ftinction  by 
slightly  rariant  graphical  means. 

(3)  To  compare  the  performance  function  as  obtained  by  the  pulse 
technique  with  the  calculated  performance  function  as  determined  from  the 
aerodynamic  deriratires  obtained  from  wind  tunnel  tests. 

(4)  To  rerecd  any  inherent  uncertainties  or  limitations  in  using  the 
pulse  technique  to  determine  performance  functions. 

To  obtain  the  data  necessary  to  inTestigate  the  abore  objectires,  the 

following  analyses  were  actually  madet 

TABU  1 

Direction  of 
Pulse  Input 


Bun  lo. 

t"^r  1 

Method 

4665 

Graphical 

4661 

f"3[«.  V 

Graphic al 

4662 

^"kt,  V 

Graphical 

4665 

C^Ci,  «,] 

Graphical 

4662 

t^'ic..  V 

Graphical 

4665 

Numerical 

'[».  w,] 


(♦)  defl. 

(-)  defl. 

(-)  defl. 

(■►)  defl, 

(♦)  defl. 

(•^)  defl. 
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In  determining  the  performance  function  graphically,  three  different 
methode  were  utilized  to  detertaine  the  amplitude  ratio  and  phase  angle  at 
each  selected  forcing  frequency.   The  first  method  utilized,  the  general 
method  as  outlined  in  Ref .  3.  that  is,  the  graphical  summation  of  input  and 
output  vectors  as  indicated  in  the  following  erpressicn: 


-jnw^T 


-Jnw-A 


^^(ln)^^n^»       ^  [FT][UT(t)]^ 


(D-1) 
n  s  1,3,3,4,5, .... 


Figure  D-1  shows  a  sample  plot  of  this  method  for  w.  «  1.25  rad/sec  using 

run  number  4665,  for  determining  the  performance  function  [^1^*)^*  v   )'      ^^ 

a  X 
this  plot ,  the  effect  of  the  lightly  damped  sinusoid  on  the  phase  angle  and 

amplitude  ratio  is  apparent.  Each  cusp  represents  one  half  cycle  of 

oscillfttlon.   In  order  for  the  amplitude  to  be  attenuated  to  0.05  of  the 

original  amplitude,  the  plot  must  "be  continued  until  a  total  of  approximntely 

6.5  cusps  have  "been  plotted.  The  majority  of  the  oscillograph  records  are 

not  long  enough  to  allow  more  than  three  or  four  cusps  to  "be  plotted,  not  is  it 

practicable  to  obtain  transient  responses  of  this  length  with  a  usable 

uncertainty  level.   Thus  the  AR  or  PA,  or  both  AR  and  PA  are  subject  to 

error  due  to  deleting  the  vector  summation  representing  several  additional 

half  cycles  of  the  oscillatory  mode.  Furthermore,  the  desired  matching  of 

the  curve  required  that  the  response  be  divided  into  0.1  sec  increments, 

requiring  that  as  high  as  84  vectors  be  plotted  (8.4  sec)  at  each  forcing 

frequency  —  a  tedious  job  resulting  in  a  rather  high  degree  of  uncertainty. 

In  order  to  simplify  the  graphical  method  by  requiring  fewer  vectors 

to  be  plotted  to  determine  the  AR  and  PA  at  each  frequency,  and  to  obviate 

the  error  Incurred  by  omitting  two  or  three  half  cycles  of  the  oscillatory 

mode,  a  second  graphical  plotting  system  was  used.   In  this  system  vectors 

were  plotted  in  Identically  the  same  manner  as  previously  discussed  until 
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the  time  respon**  was  reduced  to  a  remainder  which  was  solely  (for  practical 

purposes)  a  pure  damped  sinusoid.  At  this  point  a  single  vector  was  add»d 

to  the  response  vector  sumnation  to  represent  the  contribution  of  the  long 

"tail"  of  the  sinusoid. 

Figure  D-2  is  a  sample  plot  of  this  method  at  the  9mam   forcing  frequency 

and  run  (No.  4666)  as  used  in  Fig.  D-1.  The  accuracy  of  this  method  depends 

on  the  accuracy  with  which  the  characteristics  of  the  damped  sinusoid  may 

"be  determined.  The  characteriptlcs  were  determined  by  the  method  demonstrated 

in  Fig.  8  of  Ref.  3.  Knowing  that  the  oscillatory  motions  in  both  yawing  and 

rolling  angular  velocity  have  the  same  frequency  and  damping  ratio 

(see  eqs.  (D-IO)  and  (D-ll)  several  records  of  both  motions  were  used  to 

measure  the  period  and  peak  amplitude  ratios,  and  average  values  for  the 

parameters,  DR  and  W  ,  were  thus  obtained.  The  logarithmic  decrement  curvet 

shown  in  Fig.  7  of  Ref.  3  were  used  in  determining  DR.  The  records  agreed 

quite  closely  in  the  period  of  the  oscillation,  but  determination  of  DR 

required  the  average  of  several  records,  giving  the  following  characteristics: 

DR  «  0.16 

W  =  1.46  rads/sec 
n  ' 

In  the  case  of  the  rolling  angular  velocity  response,  it  was  found 
that  the  response  was  nearly  a  pure  sinusoid  after  the  third  peak  overshoot. 
(See  Fig.  D-6) .  As  a  consequence,  the  general  performance  function  equation 
was  rewritten,  designating  the  time  at  the  second  peak  overshoot  as  t  ,  and 


1,2, . .  »p 


the   ordinate  at   tp  as  q^.     [■  (j.  ^  r      j^  f^   o  (;rT][UT(t)] 

tn  c  i.  ,c  ,  , ,  t 

i    ., 

m  =  p,p  ▼  1 . . » 


[(£q^(T).  ^°        [yT][UT(t)]^. 

'^   ~  P«P 
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(zero  time  la^,  t  ,  is  assumed  in  the  pr«c«<iing  equation) 


The  eecond  term  in  the  numerator  may  now  be  replaced  by  r  single  rector 
representing  the  contribution  of  the  damped  sinusoid  from  t  to  infinity, 


yielding, 


2  q/   .x(t  )e 
^(out;  n 


■Jnw^i 


[rT][UT(t)] 


in  out 


'^'dnj'V* 


-jnw^ATj^ 


n  g  1 ,2 , . .  »p 


^^^^3[q<,,0(,3,)(T)]  . 


[FT][UT(t)]^ 


n  3  1,2, ..• 


(D-3) 


^^(in)^'^n^  ^ 


•Jnw^i 


^  [PT][UT(t)]^ 

n  —   l,c,o,».» 

where  At  and  At  are  the  time  intervals  used  in  the  output  and  input 
respectively. 

Any  tine  latfr  than  the  time  at  which  the  response  becomes  a  pure  sinusoid 
may  be  selected  as  the  point  at  which  the  remaining  oscillatory  motion  nay  be 
represented  by  the  Tourier  Transform  of  a  damped  sine  or  cosine  wave.  However, 
it  i«  advantageous  to  use  a  ''peak  point'*  or  a  point  where  the  sinusoid 
crosses  the  axis  in  order  that  the  phase  angle  in  the  sinusoid  expression  may 
be  eliminated.   In  these  two  cases  the  term  is  reduced  to  the  following  Fourier 


Transform  expressions: 
(t  at  peak) 


-•^0(060} 


im^ 


-DBw  t 
e      cos  w  t 

0 


3  \ 


(DR)  t\^f     |tan-\F|)-tan'^  jM^    ^ 


^j[l-   {TR)^f   *  [2(DH)(yR)]' 


e 

e 

(IK4) 


(t  at  axis) 
P 


L"]q 


-(DH)(w^)t 
oCosc)      *■   -"T)  o 


-jpw^i 


q  w 
0  o 


c 


-1 


-tan    2(DR)(rB) 


1  -  (TRff  -^   (2(DH)(rH)^ 


>  e 


■JPW^T^ 


(D-5) 
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The  time,  t   should  be  chopen  to  fall  on  one  of  the  ordlmitee  determined  by 

the  choice  of  At   to  avoid  disturbing  the  area  apDroxlnfttlon  originally 

eatabliahed  by  the  leoecelee  triangles.   If  a  peak  maxlHwin  or  mlnlmun  Ig 

selected  as  t  ,  the  area  approximation  by  triangles  is  interrupted  even 

through  t   may  lie  on  one  of  the  ordinates  of  the  fitted  triangle.   In  this 

case,  a  right  triangle  term  mupt  be  added  to  the  Fourier  Transform  of  the 

output  if  the  resulting  output  vector  is  to  be  exactly  correct.   However, 

the  initial  amplitude  of  the  cosine  wave  starting  at  t  was  small  enough 

P 

in  the  runs  analyzed  in  this  project  to  render  the  right  triangle  term 
ineffectual  in  changing  the  AR  or  PA  of  the  performance  function.  The 
method  is  demonstrated  in  Fig.  D-4. 

Figures  D-1,  D-2  and  D-3  compare  the  two  grat)hical  methods  used  in 
determining  the  performance,  function,  [^]f*w«   w  )'     Figure  I>-1  shows 
the  AH  and  PA  determined  by  plotting  the  maximum  number  of  vectors 
available  from  the  oscillograph  record.  Figure  I>-2  shows  the  same  AB 
and  PA  at  the  same  forcing  frequency  as  determined  by  representing  the 
output  divided  in  two  parts  in  the  method  Just  described.  Figure  D-3 
shows  a  comparison  of  the  AR  and  PA  results  determined  by  each  method. 
It  was  found  that  the  maximum  error  in  the  first  method,  Fig.  D-1,  appeared 
at  the  resonant  frequency. 

The  third  method  used  for  determination  of  the  performance  function 

from  pulse  response  data  is  the  method  described  in  Eef .  3,  Section  IV. 

The  oscillatory  mode  was  determined  as  previously  described,  and 

subtracted  from  the  output  response  from  t  =  0  to  o^  .  As  shown  in 

Fig.  9  of  Ref .  3,  the  remainder  of  the  output  resoonse  is  no  longer 

zero  at  t  =  T  ,  and  a  term  must  be  added  to  account  for  the  right  triangular 
o 

area  at  the  time  origin,  which  is  not  Included  in  the  area  approximated  by 
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the  lioeceles  triangles  approximating  the  remainder  of  the  response.  The 
resulting  •zpression  is  that  giren  "by  eq.  40  of  Ref.  3. 

r  vfi-i  o^OBC^ oyrem/   o     . 

n  =  1,2,.., 

-Jnw-A? 


^^''oCre.j'-.'^         [FT][(OT)(t)]^  n  =  1.2.... 

-Jnw^T 
£  [q^  (T^)]e    ^  ^  [FT][(UT)(t)]^ 


n  =  1,2,...   (D-6) 

Tlgures  3-1  and  3-2  show  the  amplitude  ratio  and  phase  angle 

respectively  of  the  perfonnance  function  for  aileron  deflection  input  and 

roll  response  output.  The  phase  angles  shown  have  "been  corrected  for  gyro 

phase  lag  as  o"btained  from  the  Fig.  1>-13.  Undamped  natural  frequencies  and 

damping  ratio  for  "both  gyros  are  found  in  Appendix  B.   In  figures  3-1  and 

3-2  three  experimental  performance  functions  are  shown  which  are  the  result 

of  analysing  graphically  three  different  oscillograph  records  (runs  4665, 

4661  and  4662)  .  The  uncertainty  in  reproducing  results  nay  be  otserved 

from  these  plots.   The  curve  shown  represents  the  average  value  of  the 

three  experimental  points  at  each  frequency.  The  greatest  deviation  from 

the  average  occurs  at  w-  =  1.25  where  the  maximum  variation  in  AH  is  5.9 

percent.  The  average  variation  is  1.6  percent.  As  noted  in  the  sample 

oscillograph  records.  Figs.  D-5  and  D-6,  no  two  input  pulses  were  identical. 

Consequently,  each  response  is  different  which  necessitated  independent 

analysis  of  each  record  .  The  reliability  of  the  pulse  input  technique  in 

reproducing  results  appears  to  he  excellent,  Judging  from  the  results  of 

this  limited  InveBtigation. 

From  the  experience  gained  in  use  of  the  analyzing  systems  previously 

outlined,  the  following  advantages  and  disadvantages  were  noted  in  applying 

these  systems  to  time  records  having  a  long  *'tail''  due  to  a  lightly  damped 

sinusoid . 
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The  only  remaining  curre  shown  in  Figs.  S-1  and  3-2  which  has  not  been 
discussed  Is  the  calculated  performance  function  uaing  wind  tunnel  stability 
derivatlres  for  the  B-25.  Due  to  differences  "between  the  model  used  in  the 
wind  tunnel  tests  and  the  actual  alrplsme,  these  derivatives  have  "been 
modified  to  more  closely  represent  the  true  derivatives  of  the  aircraft. 
The  derivatives  that  were  not  available  from  wind  tunnel  tests  were  estimated 
on  the  basis  of  previous  engineering  experience  with  comparable  aircraft 
configurations.  For  use  as  a  means  of  comparison  in  this  project,  the 
derivatives  were  corrected  to  include  the  values  of  moment  of  Inertia  about 
the  X  and  Z  axes  as  foiind  for  the  B-25  used  in  the  testing  configuration 
(See  Appendix  A) • 

The  derivatives  as  used  to  determine  a  calculated  performance  function 

are  as  follows: 

^  ^£         /^8b^     o  ,,  w 
L     =  -- — r-    —75 ■  -2.71  l/seo 

I    .  -x^-     -—-    -  ♦.673  1/.M 

K        --TT  2^  - -.00964  l/ft-.eo 

h.      =4r  HT^   = -3.175  l/..c^ 
a     ZX 

d  0       /3ovn 
H         -    -TT'        li  '  ^.00648  l/ft-seo 
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'  "-       ^VSb^ 


dC 


^(|§)     '*« 


T 


do 


f 4!^        -  --13  l/.ac 


K 


n    -^18^        «.i.66     l/.ec' 


*'  ^  *r  ^hz 


Alt.   »  10.000  ft 
Air  spd  »  175  aph 

hz 


I--     =  -1930  «ltig  ft^ 


Ijj  -  63.000  slu«  ft^ 
I22  «  120.000  slug  ft^ 
Weight  «  26.000  lb 

The  sign  oonvention  for  the  deriratires  and  equations  of  motion  and 
the  symbols  representing  displacements,  moments,  relocities.  etc..  are  in 
accordance  with  NACA  designation  system.  PositiTe  directions  along  the 
axes,  viewed  from  the  origin,  are  forward,  right  wing,  and  do%m  for  the 
X,  T.  and  Z  axes  respectively. 

The  following  assumptions  were  made  in  writing  the  equation  of  motion; 

1.  Damping  moments  due  to  control  surface  motion  are  negligible. 

2.  The  wind  sixes  and  the  aircraft  reference  axes  are  considered  to 
be  superimposed. 

3.  Ho  gyroscopic  effects  were  considered. 
Ih«  equations  of  motion  arex 

a.  Summation  of  forcts  along  T  axis 

nf  +  rfJr-  ^^▼'^^^  (^7) 

b.  Summation  of  rolling  moments 

T^.|it^^    ^p^    ikr^li^    ^    6  (D-8) 

jOC  '^     dr  d-p  or  xZ  36         a 

CI 

c.  Summation  of  yawing  moments 
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SlTldlzic  Mteh  equiition  ^7  the  appropriate  Bonent  of  inertia  or  the 
naae,  taking  the  Laplace  tranefom  and  solring  for  p  and  r,  the  follovin^ 
perforaance  functions  were  ohtainedt 

[P,](,  „  )  -  •?47P  3  (g^  *  ,S86g  *  ;,6> 

^  a  X'    (S  ♦  .377  S  •»•  1.78)  (S  ♦  2.9)  (S  -  .0013)     (B-IO) 

where  p  ■  Wj. 

[PF]        .   .184  (S  t-  3.51)(8  •»•  1.097)(8  -  .837) 

^*a  *Z^    (S^  ■••  .377  S  +  1.78)(S  ••-  2.9)(8  -  .0013)      (»-ll) 

where  r  ■  W_. 

The  calculated  perfomanoe  function  [K"]/.  w  \  ehown  In  Tigs.  3-1 

^  a  l' 
and  3-2  hare  the  eane  general  shape  as  the  experiaental  perfomance 

function.  The  BlBali^liaaat  of  the  resonant  frequencies  of  the  two  plots 

is  Indicated  in  the  difference  in  characteristics  of  the  oscillatory  modes 

of  the  calculated  and  ezperiamtal  perfornance  functions.  The  characteristics 

for  the  calculated  and  experimental  oscillatory  aodes  arei 


DB 

W 

A 

Xxperlnental 

.15 

1.45 

Calculated 

.141 

1.36 

Since  the  oscillatory  aode  Is  the  predominant  mode  in  establishing 
the  peak  ralues  of  the  perforaance  function,  the  difference  in  the 
magnitude  of  the  two  peaks  should  represent  approximately  the  same 
percentage  difference  as  indicated  hy  the  difference  in  their  damping 
ratios,  from  Tig.  3-1  the  peak  of  the  experijtental  performance  function  is 
93.8  percent  of  the  peak  ralue  of  the  calculated  performance  function,  while 
the  corresponding  ratios  of  the  two  DH  is  93  percent. 

The  calculated  performance  function  shows  a  very  small  dirergent  root. 
Howerer,  it  is  difficult  to  state  whether  the  B-25J  actuedly  is  spirally 
dirergent  or  not.  The  oscillograph  time  records  do  not  indicate  a  dirergence. 
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However,  in  order  to  deterMiae  whether  suoh  dlTergenoe  exiett,  it  would  be 
neceesary  to  take  a  record  of  at  least  15  seconds  to  detemlne  the  presence 
of  this  unstable  mode.  During  this  period  snail  perturbations  in  roll  or  yaw 
due  to  disturbances  of  the  air  nass  would  probably  nullify  the  accuracy  of 
any  such  record.  In  any  case,  for  practical  considerations,  a  root  as  snail 
as  the  calculated  performance  function  indicates,  whether  positire  or 
negative,  would  hare  negligible  effect  on  the  performance  function. 
The  sensitivities  were  conpAred  using  the  following  assumptions: 

1.  The  unstable  root  encountered  in  the  characteristic  equation  of 
the  calculated  performance  function  was  considered  negligible 
(0.0013  e  0). 

2.  The  experimental  performance  function  contains  no  divergent  roots. 
(An  assumption  implicit  in  positioning  the  base  line.) 

Using  the  first  assumption,  the  sensitivity  for  the  calculated  response  was 

determined  by  rewriting  the  i>erformance  function  as  a  sensitivity  multiplied 

by  a  frequency  function.  The  same  assumption  applied  to  the  ezperinental 

curve  gives  the  ratio  of  the  area  under  the  response  curve  in  the  tine 

domain  to  the  area  under  the  input  curve  in  the  time  domain  as  the  sensitivity 

of  the  experimental  performance  function. 

To  determine  the  area  contributed  by  the  oscillatory  mode  of  the 

response  extending  beyond  the  limits  of  the  oscillograph  record,  the  area 

was  measured  out  to  t  »  t  ,  (t  again  representing  the  point  at  which  the 

P    P 

response  beconei  a  pure  daunped  sinusoid)  and  the  curve  beyond  time  t  was 

integrated  analytically  and  added  to  the  area  from  t  =  0  to  t  .  The 

expression  thus  becomes  ^   ^^ 

*       -DW  t 

[I  q  (t  )  At  ]        ^    J  q^(T)  e    °  cos  w  t 
0  °    •  t  =^0  ° 

^(A)(6  W_)     ,    ,  ,  "  [Z  q  ?T  )  At  1 

a  X  experimental  ^1  n    l"" 

n  "^  1,2,,.. 
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The  Benaltirities  thua  determined  are: 

S(A)(fi  Wj.     .    "-^52  l/eec 
a  X  Cealc; 

^(A)(fi  W  )  *  -.945  l/tec   (Bun  lo .  4666  Input  ••♦•«) 

a  X  (exper)  * 

-  -1.080  l/sec   (Bun  No.  4665  Input  =  -*■  t   ) 

a 

«  -.985  l/eec   (Bun  Ho.  4661  Input  =  -  fi  ) 

a 

=  -v990  l/sec   (Bun  Ho.  4662  Input  =  -  6  ) 


'(A)  (6  W_),     w   X 
a  X  (experMaT) 


■-1.0  l/sec 


An  additional  sensitirity  for  comparison  with  the  abore  may  be  obtained  fro« 

Pb 
flight  teat  data  for  the  B-25,  utilizing  curres  ©^  "5fy  ▼«  *   (Bef,  4).  TroB 

these  curres,  the  sensltlTity  is, 

a  X  (rHtest  data; 

The  comparison  of  experimental  and  calculated  BensltivitieB  are  in 
close  agreement,  as  is  the  check  with  flight  test  roll  rate  data  (Bef.  4). 

The  performance  function  for  yawing  velocity  response  to  aileron  input 
is  shown  in  Figs.  3-3  and  3-4.  The  yawing  velocity  response  as  shown  in 
Fig.  D-5  appears  to  be  an  excellent  example  to  illustrate  the  liraitationi 
of  the  pulse  input  technique.  As  noted  in  Fig.  D-5,  the  response  le 
essentially  oscillatory,  the  first  lobe  (positive)  indicating  the  adverse 
yaw  effect  in  which  side  slip  has  taken  place  to  give  the  high  righting 
moment  responsible  for  the  second  lobe  (negative)  being  of  greater  amplitude 
than  the  first  lobe.  The  pattern  of  this  response  is  noted  in  detail  since 
its  effect  in  the  graphical  SLnalysis  is  important. 

The  type  of  output  response  pattern  found  in  yawing  velocity  for  aileron 
input  shows  the  same  summation  pattern  at  each  of  the  low  frequenclee .  Conr 
paring  Fig.  D-7  with  Fig.  D-2,  it  is  seen  that  the  end  point  of  each 
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gummation  is  well  away  from  the  origin,  and  that  small  errora  made  in 
plotting  will  offer  rery  little  error  to  either  the  AA  or  the  PA  of  the 
total  output  vector.  Aleo,  the  compnrieon  at  these  low  frequenciee  with 
the  AR  and  PA  of  the  calculated  performance  function  is  satisfactory  "both 
In  yaw  and  in  roll.   (See  Fig.  3-1  and  3^2  for  roll,  and  Fig.  3-3  and  3-4 
for  yaw.)   Howerer,  a«  the  frequency  increaeee,  and  Wj^t  increases  corres- 
pondingly, the  vector  summation  of  the  yaw  vectors  begin  to  "'wind  up,*' 
the  end  points  never  getting  very  far  fron  the  origin.  Note  the  difference 
in  the  sunination  as  shown  for  W«  =  1,5  rad/sec  in  Fig.  D-7  and  the  summation 
as  shown  for  W.  =  5.0  rad/sec  in  Fig.  D-8 .  This  clearly  indicates  that  any 
uncertainty  entering  the  summation,  regardless  of  its  nature,  will  impart  a 
much  higher  uncertainty  effect  to  the  PA  of  the  output  vector  (end-point  to 
origin)  than  in  the  case  where  the  end-point  is  well  removed  from  the  origin. 
The  uncertainty  of  the  results  thus  are  increased  by  two  factors  as  the 
frequency  increases, 

1.  The  large  relative  size  of  the  last  vector  in  the  sumaation  in 
comparison  with  the  resultant  output  vector,  and 

2.  The  accumulated  errors  from  measuring  ordinates  from  the  record, 
■      Consequently,  each  system  of  analysis  as  covered  in  Table  2  must  be 

analysed  and  tested  to  determine  the  relative  uncertainty  in  approximating 
the  Fourier  Transform  under  conditions  where  the  response  is  composed  in  the 

main  of  an  oscillatory  motion.  The  general  effort  is  to  establish  the  system 

k 

of  analysis  that  will  produce  the  fewest  Inaccuracies  in  approximation  and  the 


I 


I 


fewest  uncertainties  in  plotting. 

Investigation  of  the  method  listed  in  Table  2  as  (2),  the  one  tlsed 
most  extensively  in  determining  the  response  vectors  for  the  performance 
function  in  roll  for  an  aileron  input,  certain  inherent  uncertainties  are 
noted.  Again,  it  is  emphasized  that  these  uncertainties  become  increasingly 
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important  as  the  response  "becomes  more  predominately  oscillatory.   These 
uncertainties  affecting  the  magnitude  and  phase  angle  of  the  output  rector 
are  Hated  in  order  of  the  importance  of  the  error  they  present. 

1.  Uncertainty  in  determining  the  time  at  which  a  mode  has  reached  its 

peak  value  on  the  time  response,  or  the  time  nt  which  the  output  time  response 

crosses  the  "base  line.   That  is,  correct  determination  of  the  location  ef  t  * 

P 

Figures  D-5  and  D-6,  actual  oscillograph  records,  show  the  cause  for  this 
uncertainty.  The  records  are  of  too  short  duration  to  show  more  than  a  maxi- 
mum of  one  full  cycle  of  the  pure  damped  sine  wave.  Therefore,  although  the 
DR  and  W  of  the  oscillating  mode  may  he  well  eetahlished  by  comparison  of 
many  records,  the  low  flat-topped  lohee  and  the  low  slope  of  the  curve,  plus 
the  fact  that  the  noise  to  signal  level  is  becoming  quite  high  at  low 
amplitudes,  all  tend  to  make  the  selection  of  the  exact  point  at  which  the 
curve  peaks  or  crosses  the  base  line  somewhat  uncertain  within  about  +0.1 
second.  The  effect  of  such  an  error  can  be  observed  best  at  high  frequency. 
If  W  =  6  rads/sec  and  At  =0.1,  and  the  final  output  vector  is  approximately 
the  length  of  the  vector  representing  the  oscillatory  mode,  an  error  as  high 
as  30  degrees  in  PA  may  exist.  As  the  end  point  moves  further  from  the 
origin,  this  possible  maximum  error  reduces  accordingly. 

2.  Uncertainty  in  determining  DH  and  W  of  the  oscillatory  mode. 

3.  Uncertainty  in  determining  the  amplitude  of  the  oscillatory  mode. 
On  the  basis  of  the  above  uncertainty  at  high  frequency  in  adding  a 

vector  representing  the  Fourier  Transform  of  the  oscillatory  mode,  it  may 

be  concluded  that  the  greater  the  time,  t  ,  the  less  the  uncertainty.   This 

may  be  observed  in  eq .  (D-3) .   To  simplify  that  expression  consider 

At  =  At.  ,  then 
o     i 
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£  q      (t    )•  •     n  =  l,2,...p 

[FT]  «    ft S- 


(    )(q.q^)  -Jnw  At 

^  ""        £<iMt  ).       ^    ^ 

1       n 


&  =  1  ,c  ,3 1  . . . 


1       "JP^'^-W^ 


[FT][UT(t)]^         S  qi(T^)e         ^ 


X  I  lO  g  fl  •  • 


The  larger  t  'become«,  the  greater  the  percentage  of  the  graphical  emalysii 
P 

that  Is  accompllehed  with  the  same  degree  of  approximation  in  both  nunerator 
and  denominator.  If  the  oscillograph  record  were  long  enoxigh,  the  test  results 
at  high  frequency  would  be  obtained  by  continuing  summation  of  the  output 
vectors  until  such  time  as  the  end  point  is  clearly  indicated.  This  allows 
the  same  degree  of  approximation  in  obtaining  both  the  outpcit  rector  and 
the  inuut  vector  throughout  the  entire  summation.  However,  in  addition  to 
the  tedium  of  plotting  a  large  number  of  vectors,  uncertainties  are  intro- 
duced by  the  oscillograph  record  in  picking  off  ordinates  well  out  in  the 
time  response.  The  noise  to  signal  ratio  becomes  high,  and  drift  of  the 
baseline  appears  to  be  one  of  the  evils  of  the  long  record. 

The  most  uncertain  of  the  graidiical  methods  for  use  with  the  highly 
oscillating  type  response  is  the  method  which  subtracts  the  pure 
oscillatory  mode  from  the  entire  response.  In  this  case,  the  subtracted 
oscillatory  mode  represents  a  high  percentage  of  the  total  response.  This 
large  portion  of  the  response  is  represented  by  a  single  vector  which  is  an 
exact  Courier  Transform.  Consequently,  the  output,  in  the  extreme  case  where 
the  response  is  almost  entirely  oscillatory,  is  represented  exactly,  and  the 
input  Is  approximated.   The  result  is  that  at  high  frequency,  the  output 
reaches  a  certain  fixed  phase  angle,  due  to  the  predominance  of  the  oscillatory 
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Bod«,  and  the  input  continues  to  shift  phase  angle  rapidly.  Therefore, 

it  appears  that  at  rery  high  frequency  the  oscillatory  part  of  the  output 

will  erentually  lead  the  input.  The  lagging  phase  of  the  remainder  should 

theoretically  compensate  in  such  a  manner  that  the  vector  representing  the 

total  performance  function  will  hare  the  proper  phase  angle.  But,  the 

rector  representing  the  remainder  is  rery  small,  and  therefore,  although 

its  phase  lag  nay  "be  increasing  In  a  like  manner  to  the  denominator  (input), 

its  contribution  to  the  phase  angle  of  the  vector  summation  is  very  small. 

Furthermore,  the  additional  manipulation  of  the  time  response  data  in 

obtaining  a  "breakriovn  into  components  offers  added  uncertainty  in  the  vector 

summation. 

Sample  records  are  provided  to  show  the  vector  summation  using  method 

(3)  in  which  the  oscillatory  mode  is  removed  from  t  =  0  to  t  =  00  and 

represented  by  an  exact  Tourier  Transform,  and  method  (2)  in  which  the 

Fourier  Transform  of  a  sine  wave  starting  at  t  is  added  to  the  output 

summation  for  t  »  0  to  t  =  t  .  Figures  D-9  and  D-10  respectively  show  these 

P 

two  cases.  Note  that  in  Fig.  D-9  the  rector  representing  the  oscillatory 
mode  is  the  largest  of  the  three  vectors  comprising  the  total  response 
vector.  Furthermore,  it  has  attained  its  maximum  lagging  phase  angle  to 
within  about  5  degrees.  The  phase  angle  from  the  performance  function  at 
W.  =  6  by  this  method  is  -49  degrees  (corrected  for  gyro  phase  amgle 
from  Fig.  D-14) .  Figtire  D-10  shows  the  results  of  using  the  same  approxi- 
mation in  both  input  and  output  over  the  greater  part  of  the  response,  method 
(2) .  The  addition  of  the  vector  representing  the  sine  wave  summation  is 
nearly  negligible  in  its  effect  on  the  total  response  vector.  The  phase 
angle  at  this  same  frequency  by  this  method,  (2),  is  -87.5  degrees. 

As  previously  discussed,  adding  the  sine  term  involves  uncertainty  in 
properly  locating  t  ,  This  error  may  be  reduced  to  a  negligible  size  or 
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eliralnnted  "by  the   following  meguis.      Select   t     where   the  output   croBses  the 

P 

base  line  and  ohserve  the  magnitude  of  the  vector  resulting  in  comparison 

with  the  magnitude  of  the  last  vector  in  the  summation  from  t  =  0  to  t  =  t  . 

P 

If  the  vectors  are  approxinBtely  the  same  size,  no  appreciable  error  will 

result  from  an  error  in  locating  t  .   If  the  vector  representing  the  added 

sine  term  is  large  compered  to  the  last  vector  of  the  summation,  t  should 

be  moved  to  the  next  null  point  of  the  time  response,  which  will  decrease 

-W  (DR)(T/2) 
the  magnitude  of  the  sine  wave  amplitude  by  the  factor  e  • 

It  is  the  opinion  of  the  authors  that  the  added  sine  wave  method  of 

component  summation  is  subject  to  the  least  inaccuracies,  of  the  methods 

investigated,  in  determining  the  performance  function  of  a  response  which  is 

essentially  oscillatory.  The  method  affords  the  sane  degree  of  approximation 

in  both  the  input  and  response,  until  such  time  as  the  response  is  subjected 

to  inaccuracies  of  noise  and  instrumentation.  At  this  point  (t  )  the  small 

P 

remainder  is  represented  by  an  exact  transform. 

Referring  again  to  Figs.  3-3  and  3-4,  the  effects  of  various  analysis 
methods  in  obtaining  the  performance  function  are  sho%ra.  The  method  numbers 
are  those  listed  in  Table  2.  As  previously  noted,  the  variation  in  PA  from 
the  calculated  performance  function  at  high  frequencies  is  very  high,  and  it 
is  difficult  to  explain  the  variation  in  terms  of  the  dyneonics  of  the  air- 
plane. The  assumptions  listed  in  setting  up  the  theoretical  equations  of 
motion  for  the  test  vehicle,  limit  the  phase  shift  to  90  degrees  of  lag  at 
high  frequencies.  However,  refinements  which  were  excluded  by  these 
assumptions  aay  add  about  10  to  15  percent  additional  phase  lag  at  high 
frequency,  as  shown  in  the  theoretical  performance  function  found  in  Ref .  7. 
The  performance  function  applies  to  the  B-26,  an  aircraft  similar  in  dimensions 
and  mass  configuration  to  the  B-25J .   The  experimental  PA,  at  high  frequencies 
are  in  disagreement  with  the  calculated  performance  functions  determined  both 

for  the  B-25  and  the  B-26.  However,  the  very  low  AH  at  the  frequencies  of 
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uncertain  phase  angles  makes  PA  uncertainty  relatirely  \inlmportant  In 
deterainlng  the  performance  functions. 

The  scattering  of  points  In  the  doubtful  phase  angle  range  shows  the 
effect  of  various  methods  of  determining  the  transform  of  the  response.   It 
was  found  "by  halving  the  At  ,  the  phase  angle  was  negligibly  affected. 
Compare  Fig.  D-8  with  Fig.  I>-11.  This  indicates  that  V^t  has  not  reached 
a  critical  value.  (Ref.3). 

The  eensitivities  for  this  performance  function  were  determined  in  the 
same  manner  as  used  in  the  performance  function  for  rolling  velocity  response 
to  input  aileron  responses. 

As  previously  noted,  the  calculated  performance  function,  (eq.  D-ll) , 
has  a  very  small  root  in  the  denominator.  Such  a  small  root  either  positive 
ot   negative  would  not  be  apparent  in  the  time  response  and  would  be  Impllcltiy 
considered  to  be  eero  by  the  action  of  drawing  In  the  base  line.  To  be 
consistent  for  comparative  purposes,  the  same  assumption  ntist  be  made  with 
the  calculated  performance  function. 

Thus,  to  get  a  valid  comparison  let  Jw  -  0.0013  «  Jw  in  the  calculated 

performance  function.  Then  by  writing  the  performance  function  as  a 

sensitivity  multiplied  by  a  frequency  function,  a  sensitivity  that  can 

be  compared  with  the  ratio  of  the  areas  under  the  experimental  output  and 

input  curves  is  obtained.  The  comparison  is: 

S/A\r«  w  T     ■  0.209  l/sec 
^^^•^^  ^Z^calc 

^^aUa  w  1      "  ''^^'^   l/sec      (from  run  no.  4665  •♦•  fi   input) 
U;L6^  *Z-»exper  * 

S/.\r«  «  T      ■  •1'''2  l/sec      (from  run  no.  4662  -fi   input) 
VAno  w_l  a 

•^  a  Z-'erper 

In  general,  the  results  obtained  bring  out  several  points  concerning 

the  application  of  the  pulse  technique  as  a  means  of  determining  the 
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performance  functions  of  dynamical  eyetems.  B«8ult«  obtained  in  rolling 
Telocity  response  for  aileron  inputs  showed  very  close  agreement  with  the 
theoretical  performance  function  throughout  the  entire  range  of  frequencies 
investigated.  The  comparison  of  time  response  patterns  or  conformations 
indicate  the  response  pattern  involving  the  least  inaccuracies  in  analysis 
is  one  which  responds  quickly,  reaahing  its  maximum  amplitude  in  the  early 
part  of  the  time  response.  Under  these  conditions,  the  end  point  will  tiaiiallx 
■be  well  removed  from  the  origin,  and  uncertainties  in  added  vectors  represent- 
ing the  Tourier  Transforms  of  oscillatory  modes,  right  triangles,  etc.  will 
have  a  minimum  effect  on  the  AH  and  PA  of  the  performance  function,  until 
very  high  frequencies  and  small  AR  prevail. 

▲  potential  source  of  error  lies  in  establishing  the  true  baseline  for 
each  response  and  input  time  record.  An  error  In  the  location  of  this  line 
produces  a  double  error  in  the  results  in  that  as  much  false  area  is  added 
to  one  side  of  the  line  as  is  subtracted  from  the  other  side  of  the  line.  In 
the  case  of  records  having  the  long  tail  characteristic  of  a  lightly  damped 
oscillatory  mode,  the  envelope  should  be  drawn  in  determining  the  base  line, 
particularly  if  there  is  any  drift  in  the  center  of  the  response  record.  This 
is  demonstrated  in  Fig,  8,  Ref.  3.   In  any  case,  measurements  from  the  reference 
trace  to  the  center-line  should  be  made  at  short  intervals  and  connected, 
rather  than  passing  a  line  between  two  measured  points  at  the  extremities  of 
the  record.  This  corrects  for  warping  of  the  record  incurred  in  the  drying 
process.   If  nominal  care  is  exercised  in  laying  out  the  baseline,  and  in 
picking  off  ordinates  from  the  record,  the  uncertainties  in  reproducing 
results  will  be  satisfactorily  low. 
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FLIGHT  TEST  HECOBDS 
FLIGHT  NO.  1 
Date:   28  July  1949  Duration:  2.5  hri 

Airplane:  B-25J  Ho.  44-30328 
Crew: 

Pilot:  C.  0.  BoBtroM 
Co-pilot:  F.  H.  Michaelie 


i 


ObBerrere:  J.  C.  Woetton 

A.  A.  Hollander 

Purpose  of  flight: 

Instrument  shakedovm 
Airplane  configuration: 

Gross  Weight  -  26,100 

C.  G.  Percent  MAC  -  27.4 
Instrumentation  installed  for  recording  elevator  position,  rudder  position, 
rate  of  roll  and  rate  of  yaw. 
Flight  Procedure: 

(1)  Trimmed  aircraft  at  10,000  ft  PA  and  175  mph  IAS. 

(2)  Ohserved  action  of  oscillograph  of  input  and  response  traces 
using  test  pulses  of  positive  and  negative  aileron  and  rudder  displacements. 

(3)  Obtained  records  of  positive  and  negative  ailerons  and 
rudder  pulses . 

Results  and  Coonents: 

Frequency  regulation  of  the  inverter  was  very  poor.  Frequency  varied 
from  380  cps  to  405  cps .  Noise  level  on  the  oscillograph  records  was  so 
high  that  the  records  were  not  satisfactory  for  reduction. 
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FLIGHT  NO.  2 
Date:   2  August  1949  Duration:   1.5  hr« 

Airplane:  B-25J  No.  44-30328 
Crew: 

Pilot:  C.  Collins 

Co-pilot:  J.  B.  Bain 

Oteerrort*  T.  H.  Michaelis 
J.  C.  Wootton 

Purpose  of  Flight: 

To  o'btain  oscillograph  records  of  airplane  response  to  aileron  and 
rudder  inputs. 
Airplane  configuration: 

Gross  weight  -  26,100 

C.  G.  Percent  MAC  -  27  .4 
Flight  Procedure: 

(1)  Trimmed  at  10,000  ft  PA  and  175  mph  IAS. 

(2)  Obserred  action  of  oscillograph  treicet.for  test  pulses  of  posltire 
and  negatiTe  aileron  pulses. 

(3)  Obtain  records  of  positire  and  negative  aileron  and  rudder  pulses. 
Results  emd  Connents: 

No  records  obtained  due  to  aalfunctioning  of  recording  equipment. 
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FLIGHT  HO,   3 
Date:      2  August   1949  Duration:      1.5  hra 

Airplane:      B-25J  No.  44-30328 
Crew: 

Pilot:     C.  Collins 

Co-pilot:     J.  B.  Bain 

Obserrers:  F.  H.  Mlchaells 
J.  C  .  Wootton 

Purpose  of  Flight: 

To  obtain  oscillograph  records  of  aircraft  response  to  aileron  and 
rudder  Inputs. 
Airplane  configuration: 

Gross  weight  -  26,100 

CO.  Percent  MAC  -  27 .4 
Procedure  this  Flight: 

(1)  Trlmned  at  10,000  ft  PA  and  175  mph  IAS. 

(2)  Ohsenred  action  of  oscillograph  traces  for  test  pulses  of  posltlre 
and  negative  aileron  pulses. 

(3)  Obtain  records  of  positive  and  negative  aileron  and  rudder  pulses. 
Results  and  Comaents: 

The  records  obtained  had  an  unacceptable  noise  level,  but  were  somewhat 
improved  over  preceding  flights.   Input  pulses,  while  acceptable,  were  too 
long  for  ease  of  data  reduction. 
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7LiaHT  10.  4 
Date:     5  August  1949  Duration:      1.50  hr 

Airplane:     B-25J  Ho.  44-30328 
Crow: 

Pilot:     Capt.  Wanrlck 

Co-pilot:     J.  B.  Bain 

Obsenrers:  J.  C.  Wootton 

r.  H.  Miehaelii 

Purpose  of  Flight: 

To  o'btain  oscillograph  records  of  airplane  response  to  aileron  emd 
rudder  displacements* 
Airplane  configuration: 

Gross  weight  -  26,100 

CO.  Percent  MAC  -  27 .4 
Replaced  faulty  inrerter,  and  installed  new  Lord  shock  mounts  on  gyro  can, 
Tllght  Procedure: 

(1)  Trimmed  at  10,000  ft  PA  and  175  mph  IAS. 

(2)  Obserred  action  of  oecillograph  traces  for  test  pulses  of 
posltlre  and  negative  aileron  pulses. 

(3)  Obtain  records  of  positire  and  negatire  aileron  and  rudder  pulses. 
Besults  and  Comnenta: 

Flight  discontinued  due  to  wiring  difficulties. 
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7LIGHT  NO.  5 
Date:   5  August  1949  Duration:   1.5  hre 

Airplane:   B-25J  No.  44-30328 
Crew: 

Pilot:  Capt .  Warwick 

Co-pilot:  J.  B.  Bain 

Obeerrers:  J.  C.  Vootton 

T.   H.  Mlchaells 

Airplane  configuration: 

Grose  weight  -  26,100 

C.  G.  Percent  MAC  -  27 .4 
Changes  since  last  flight: 

Loose  connection  in  300  volt  system  found  and  repaired. 
Flight  Procedure: 

(1)  Trimmed  at  10,000  ft  PA  and  175  mph  IAS. 

(2)  O^serred  action  of  oscillograph  traces  for  test  pulses  of 
positive  and  negatire  aileron  pulses. 

(3)  Obtain  records  of  positive  and  negative  aileron  and  rudder  pulses 
Results  and  Comments: 

The  records  obtained  were  considered  to  be  useable.   The  noise  was  of 
insignificant  amplitude,  and  the  input  pulses  were  of  the  correct  size, 
however,  the  roll  gyro  picked  up  a  vibration,  apparently  from  its  mount, 
signal  of  about  10  cps  which  would  make  data  reduction  somewhat  difficult. 
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FLIGHT  NO.  6 
Date:   9  Axiguet  1949  I>aratlon:   1.5  hre 

Airplane;   B-25J  No.  44-30328 
Crew: 

Pilot:  C.  Collins 

Co-pilot;  J.  B.  Bain 

Observers;   F.  H.  Michaelis 
7.   Smith 

Airplane  Configuration: 

Grose  weight  -  26,100 

CO.  Percent  MAC  -  27.4 
Changes  since  last  flight: 

Intalled  stop  to  remove  shock  mounts  while  airborne  and  wedged  gyro 
container  rigidly  to  airplane  to  tllminat*  rihration  of  the  mount  from 
being  recorded. 
Prootdure  this  flight: 

(1)  Trimmed  at  10,000  ft  PA  ajid  175  mph  IAS. 

(2)  Observed  action  of  oscillograph  traces  for  test  pulses  of  positive 
and  negative  aileron  pulses. 

(3)  Obtain  records  of  positive  and  negative  aileron  and  rudder  pulses. 
Results  and  Comoents: 

The  records  obtained  on  this  flight  were  very  good.   The  attenuator 
settings  used  were: 

Function  Attenuator 

Wj  4 

6  5 

a 

fi  e 

r 
These  attenuations  allowed  the  traces  to  remain  on  the  pa|)er  and  gave  large 

enough  traces  to  be  readable,  light  records  were  taken  and  static 
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calibrations  with  the  master  pickoff  were  made  in  flight. 
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YAWING    PATE 


'•rm^^  SJ^  t^    *iii"'<^>^i        ^^ 
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APPINDII  I 
OLOSSABT 
The  eyeten  of  notation  adopted  for  use  in  this  thesis  has  been 
selected  l)eoHU9S  of  its  peculiar  usefulness  in  representing  opernting 
system  performance.   In  particular,  it  has  the  following  specific 
advantages: 

1.  It  is  easily  learned. 

2.  It  is  adopted  to  a  wide  range  of  situations. 

3.  It  is  "built  up  almost  exclusively  of  characters  found  on  the 
keyboard  of  a  standard  American  typewriter. 

4.  Any  one  of  the  compound  symbols  of  the  system  is  readily  inter- 
preted without  recourse  to  an  extensive  glossary. 

A  representative  list  of  the  primary  symbols  and  those  used  as  sub- 
scripts is  given  in  the  following  tabl«. 
Primary  Symbolf 

A  Angl« 

(ah)  y  x/        X  Input  -  output  amplitude  ratio  for  a 


in  out  given  operating  component  =  q/  ^\ 


a 


^(in)a 
c  Linear  viscous  damping  coefficient 

(CT)  Characteristic  time 

(DR)         .  Damping  ratio  =  — ° 

2  VS/   N  m 
\   (er; 

f 
(FH)  Frequency  ratio  =  tt- 

n 
(FT)  Fourier  transform 

g  Force  of  gravity/unit  of  mass 


I-j  Moment  of  inertia  of  airplane  about 

X  axis 

Ij2  Product  of  inertia 
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^zz 

L 

L 
P 

L 

r 

h. 

L 

T 

(LT)i:q(t 

)] 

Moment  of  Inertia  of  airplane  about 
Z  axis 

Boiling  moment 

Rolling  moment  due  to  rolling  velocitj 

Boiling  moment  due  to  yawing  velocity 

Boiling  moment  due  to  aileron  displacement 

Boiling  moment  due  to  side  slip  Telocity 

Direct  Laplace  transform  of  q(t)  »  (i(S) 

■  Mass 

H  Tawing  moment 

N  Tawing  moment  due  to  rolling  velocity 

P 

N  Tawing  moment  due  to  yawing  velocity 

r 

H  Tawing  moment  due  to  side  slip  velocity 

N,  Tawing  moment  due  to  aileron  deflection 

oa 

Ng  Tawing  moment  due  to  rudder  deflection 

^P^/  w        V  Input  -  output  performance  operator  for 

^In  ^out  a  given  operating  component 

(Pa) /  w        V  Input  -  output  dynamic  phase  angle  for  a 

*^in  out  given  operating  component 

(PF)/  >.  /        X  Input  -  output  performance  function  for  a 

in  out  given  operating  coarponent,  where 

1(ln)  -  "(in).  • 

^(out)  '  ^(out)a  * 

q/   \  Input  quantity 

q,   V  Output  quantity 

^Cout; 

q(t)  A  quantity  which  is  a  function  of  the 

real  variable,  time 

BJ  Bight  triangle 

fl  Complex  variable  used  in  Laplace  transform 

theory 

Static  sensitivity  for  a  given  operating 


'(  )U.«  q^„J 


in  ^out  component 
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t^  Time  delay 

At  Finite  increment  of  time 

U  Airplane  forward  Telocity 

(OT)  Unit  triangle 

(tJBf)  Unit  right  triangle 

T  Acceleration  along  T  axle 

W/  \  Airplane  angular  velocity  about  axia  (  ) 

W  Angular  forcing  frequency 

W  Natural  frequency 

T  Side  force 

T  Side  force  due  to  side  slip  velocity 
A  dot  over  a  variable  represents  its  derivative  with  respect  to  time 

(..g.  (■).  -4^). 

Two  dots  over  a  variable  represents  a  second  derivative  with  respect  to 
time  (e.g.  (   )  «  - — ^  ). 

The  followlTig  Oreelc  letter  sjnnbols  are  retained  because  their  meanings  have 

become  thoroughly  identified  with  these  symbols  in  the  aeronautical  engineering 

field. 

6  Total  aileron  deflection  (sura  of  left  and 

right  aileron  deflection)  (positive  deflection 
of  right  aileron  is  downward) 

t  Rudder  deflection  (left  rudder  deflection  is 

positive) 

0  Angle  of  bank 

3  Angle  of  yaw 

tt)  Angular  forcing  frequency 
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Modifying;  symbols 

(sub: 

scripts) 

a 

Aileron 

(•) 

Amplitude 

A 

Airplane 

(dlT) 

Dirergent 

(in) 

Input 

(ind  878 ) 

Indicating  system 

i" 

Measuring 

0 

Output 

(out) 

Output 

(osc) 

Oscillating 

p 

Bate  of  change  of  bank  angle  ^  W- 

(rem) 

Remainder 

r 

Rat©  of  change  of  heading  =  V 

T 

Side  slip  velocity 
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